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Abstract: The "expansion by regions" is a method of asymptotic expansion developed 
by Beneke and Smirnov in 1997. It expands the integrand according to the scahng pre- 
scriptions of a set of regions and integrates all expanded terms over the whole integration 
domain. This method has been applied successfully to many complicated loop integrals, 
but a general proof for its correctness has still been missing. 

This paper shows how the expansion by regions manages to reproduce the exact result 
correctly in an expanded form and clarifies the conditions on the choice and completeness 
of the considered regions. A generalized expression for the full result is presented that 
involves additional overlap contributions. These extra pieces normally yield scaleless inte- 
grals which are consistently set to zero, but they may be needed depending on the choice 
of the regularization scheme. 

While the main proofs and formulae are presented in a general and concise form, 
a large portion of the paper is filled with simple, pedagogical one-loop examples which 
illustrate the peculiarities of the expansion by regions, explain its application and show 
how to evaluate contributions within this method. 
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1 Introduction 

When loop integrals involve many different scales from masses and kinematical parameters, 
it can be hard or even impossible to evaluate them exactly. The integrand may be simplified 
before integration by exploiting hierarchies of parameters and expanding in powers of small 
parameter ratios. When these expansions are done naively, neglecting their breakdown in 
certain parts of the integration domain, new singularities may be generated and important 
contributions to the full result can be missed. A proper treatment requires sophisticated 
methods of asymptotic expansions. One of them is the so-called "strategy of regions" or 
"expansion by regions" developed by M. Beneke and V.A. Smirnov [1]. The recipe for 
applying the expansion by regions to a loop integral reads as follows [1-4]: 

1. Divide the space of the loop momenta into various regions and, in every region, 
expand the integrand into a Taylor series with respect to the parameters that are 
considered small there. 

2. Integrate the integrand, expanded in the appropriate way in every region, over the 
whole integration domain of the loop momenta. 

3. Set to zero any scaleless integral. 

The sum of these contributions yields the full result of the loop integral in an expanded 
form. This recipe will be illustrated in the examples of the following sections. 

However, despite the successful application of the expansion by regions in many loop 
calculations, each of the steps in the recipe above raises questions: 

• In step 2 all expanded terms are integrated over the whole integration domain, ne- 
glecting the domains of convergence of the expansions. This generally leads to new 
singularities which obviously must be cancelled between the contributions of the var- 
ious regions. How is this cancellation ensured? 

• In the original integral each point of the integration domain contributes exactly once. 
According to step 2 each point of the integration domain contributes once per region. 
How can this double- or multiple-counting of contributions be correct? 

• How do we have to choose the regions in step 1? And how do we know that the 
chosen set of regions is complete? 

• Although it seems natural to eliminate scaleless integrals when using dimensional 
or analytic regularization, we can ask: What is the role of the scaleless integrals in 
step 3? 
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These questions will be addressed in the current paper. 

The developers of the expansion by regions have introduced their method using exam- 
ples and "heuristic motivations" (e.g. related to analogies between the regions and degrees 
of freedom in effective theories). They have shown the validity of the method for some types 
of expansions (in particular Euclidean- type limits like off-shell large-momentum expansion 
or large-mass expansion) through the agreement with existing and proven expansions by 
subgraphs [4]. They admit, however, that they cannot give a general mathematical proof 
of their prescriptions [1] and that "it is not guaranteed that expansion by regions works in 
all situations" [4]. 

A practical problem in the application of the expansion by regions is the correct choice 
of the regions. In the original paper [1] the authors determine the relevant regions from 
the structure of the poles of the propagators in the loops. They close the contour of the 
integration over the zero-component of the loop momentum and study its scaling at the 
residues depending on the size of the spatial components. In general, relevant regions can 
often be found by looking at the structure of the integrand and at singularities which arise 
in the given parameter limit. It does not matter to consider more regions than necessary: 
The irrelevant regions will simply produce scaleless contributions which are set to zero. 
The tricky point here is to avoid double-counting of regions which look different but yield 
equivalent expansions. 

Alternatively, the expansion by regions has been applied to the alpha-parameter rep- 
resentation of loop integrals [3, 4]. The double-counting of equivalent regions should not 
occur here, but there exist at least some regions (in particular the "potential region" in 
threshold expansions) which cannot easily be identified in the language of the alpha pa- 
rameters. Based on this approach, an algorithm for finding the regions in alpha-parameter 
representations has recently been developed [5]. In the present paper I stick exclusively to 
the version of the expansion by regions which is applied directly to loop integrals because 
I find this original variant of the method more natural. However, the formalism which 
I present is general and can be applied to the asymptotic expansion of many types of 
integrals, including alpha-parameter representations. 

An important check for the completeness of the regions is the cancellation of singular- 
ities. The sum of all contributions must not be more singular than the original integral, 
and if additional singularities persist, then probably the contribution from another region 
is missing. This check is, however, not sufficient to guarantee the completeness of regions, 
because some regions or a subset of them can also yield non-singular contributions such 
that their absence could remain undetected. 

Another possibility of finding all relevant regions for an integral in a given parameter 
limit is the one which I have successfully used in many calculations (e.g. [6-9]): Evaluate 
the full integral with the propagators raised to generic powers in terms of as many Mellin- 
Barnes representations [10] as needed. Extract the leading-order asymptotic expansion by 
closing those Mellin-Barnes integrals which involve the expansion parameter and taking 
the residue of the first pole next to the contour. This yields a sum of terms reproducing the 
exact result up to higher powers of the expansion parameter than those already present. 
The terms often still contain Mellin-Barnes integrals and may be too complicated to eval- 
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uate them further. But each term is characterized by a homogeneous dependence on the 
expansion parameter raised to some power which is a function of the space-time dimen- 
sion d = 4 — 2e and of the propagator powers. An examination of each term's dependence 
on the expansion parameter is usually able to tell the corresponding region needed for pro- 
ducing this contribution. If the asymptotic expansion by Mellin-Barnes representations 
is performed correctly, the prescription described here yields all the relevant regions for a 
subsequent application of the expansion by regions. 

The correspondence between the contributions of the expansion by regions and the 
asymptotic expansion via Mellin-Barnes representations is, of course, as heuristic as the 
derivations and justifications known so far for the expansion by regions. In the present 
paper I follow a different approach: I start from a general integral for which an expansion 
is required and transform the expression step by step, in a mathematically well-defined 
way, into a sum of expansions which can be identified with the ones originating from the 
expansion by regions. The resulting expression contains additional overlap contributions 
which are absent in the usual prescription for the expansion by regions. While these extra 
pieces normally yield scaleless integrals which are consistently set to zero according to 
step 3 above, there are cases — depending on the choice of the regularization scheme — 
where these overlap contributions are present and required for the correct result. Such 
overlap contributions have also been introduced in effective-theory treatments as "zero-bin 
subtractions" (see e.g. [11, 12]). 

The basic idea of the formalism which I develop in this paper goes back to a one- 
dimensional toy example [13] (see also section 3.2 of [4]). I do not treat this one-dimensional 
toy integral, but start directly with a simple d-dimensional loop integral in section 2 before 
presenting a first version of the formalism for general integrals in section 3. A generalized 
version of the formalism is elaborated in section 5, where one restriction of the first version 
is relaxed. Examples which illustrate the application of the formalism and explain the eval- 
uation of loop integrals with the expansion by regions are shown in sections 4, 6 and 7. The 
last of these examples demonstrates the relevance of overlap contributions under certain 
circumstances. The main statements and formulae are summarized and discussed in the 
conclusions of section 8. Details from the evaluation of the expanded loop integrals have 
been shifted to appendix A. Finally appendix B treats an example with a finite integration 
boundary which illustrates a subtlety raised in section 3. 

The order of the general parts (sections 3, 5 and 8) and illustrating examples (sections 
2, 4, 6 and 7) has been chosen from a pedagogic viewpoint in order to facilitate the un- 
derstanding especially for readers who are not very familiar with the expansion by regions. 
More experienced readers may skip the examples and concentrate on the general sections 
for a quick study of the main statements. The examples, however, also show how the 
general formalism can actually be applied to loop integrals and how its conditions on the 
regions are checked. Later examples use notations and conventions introduced in earlier 
examples. 
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Figure 1. Two-point loop integral for off-shell large- momentum expansion. 

2 Example: off-shell large-momentum expansion 

The two-point one-loop integral of the first example is defined by the expression 

F = ^ DA;/ (2.1) 
with the integrand / = /1/2, the propagators 

^1 = ((fc + p)2)"i = (A;2 + 2A: • p + p2)ni ^2 = _ j„2)n2 (2-2) 



and the integration measure 



Bk ^ ^^'^ e^''- I ^ , (2.3) 



where d = 4 — 2e is the space-time dimension, /i is the scale of dimensional regularization, 
and 7e ~ 0.577216 is Euler's constant. The usual infinitesimal imaginary part in the 
Feynman propagators is understood and has been dropped in the notation for brevity. 

The integral F depends on the external momentum p and the mass m. We are inter- 
ested in the off-shell large-momentum limit 

and look for an asymptotic expansion in powers of ir? jp^ . The integral also depends on the 
propagator powers n\ and 77-2, and we focus particularly on the case ni = 1, 722 = 2 which 
is depicted in figure 1 and for which the integral is finite for both small and large k. Exact 
results for this integral can easily be obtained and expanded in m? /p^ for comparison with 
the asymptotic expansion. 

In the following two subsections we first treat this loop integral according to the recipe 
formulated in section 1 for the expansion by regions before proving the correctness of this 
approach by independent mathematical transformations of the integral. 

2.1 Expansion by regions in the large-momentum limit 

The first propagator in (2.2) is characterized by the large momentum p, whereas the second 
propagator is characterized by the small mass m. It is therefore natural to assume that 
the two regions of relevance to this problem are 
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• the hard region {h), where k ^ p, 

• and the soft region (s), where k ^ m. 

By A; ~ p we mean that ah components of the loop momentum k are of the order of ( "scale 
like") the (overall size of) the momentum p, and similarly for k ^ m. Instead of dividing 
the integration domain into subdomains with explicit boundaries, the regions simply define 
scaling prescriptions for the loop momentum on the basis of which we are able to perform 
expansions of the integrand. 

For the hard region we consider j/c^l ^ and expand the second propagator in (2.2), 
while the first propagator remains unchanged: 

/.^TC)/../,, /,^T('")/..2T,<''>/. = f (2.4) 

where T^^'^ is the expansion operator of the hard region, and Tj'^^ generates its j'-th order 
expansion term. For the soft region we consider \k'^\ <^ \p'^\ and |2A; -pj <^ permitting 
the expansion of the first propagator, while the second remains unchanged: 

T ^rWr -\^T^'h -\^t(^) t - V" K)ii2 {-k'^y^ {-2k ■ py^ 

h ^ >h = 2^r. /i = 2^ T.^j.A - 2^ — y-|- — ( 2)ni+ii2 — ' 

j ji,j2 ii,i2=o 

h ^ T^'^h = h , (2.5) 



where 

ja(3- = ja+ 3p + ■■■ (2.6) 

is) 

is introduced as a shorthand notation, which is also used for other symbols. Here 
generates the j-th order expansion term of the soft region, and it is a priori not clear how 
j relates to ji and j2 because k"^ and 2k-p involve different powers of the soft momentum k. 
So we postpone this question until after the loop integration when we will rearrange the 
summation over ji and j2. Both the hard and the soft expansions have employed the 
generic Taylor expansion 



if < |X|Vi, 

(2.7) 

with the Pochhammer symbol (a)j = r(a + j)/r(a). 

According to step 2 of the recipe in section 1, each expanded term has to be integrated 
over the whole integration domain. The j-th order contribution of the hard region reads 

Ff ) = [Bk rf )/ = (my [ : , (2.8) 

J J J ' J ((A;+p)2)"^(A:2)"2+j ' ^ ' 



and the soft-region contribution with indices ji , j2 is given by 

pis) _ [^^ ^is) J _ (-ly- fj,, {kyH2k-py^ 
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These new integrals are simpler than the original integral (2.1): The hard contribu- 
tions (2.8) are massless integrals and functions only of p^. And the soft contributions (2.9), 
once each scalar product in the numerator has been separated into k ■ p = p^k^, are mas- 
sive tadpole tensor integrals of rank j2 and functions only of m^. These are all one-scale 
integrals, and it is already clear at this step from a dimensional analysis that the hard and 
soft contributions are homogeneous functions of rr? j]? with 

Ff^ OC {m^y f-"i2-J- , fW^.^ oC (m2)f-2+ii+i2/2 |^2|-ni-ii-,V2 _ (2.10) 

We also see that the expansions generate new singularities: The hard contributions (2.8) 
will become infrared-singular (for — t- 0) due to the increasing number of /c^-terms in the 
denominator, and the soft contributions (2.9) will become ultraviolet-singular (for k — )• oo) 
with more and more powers of the loop momentum in the numerator. In the particular 
case ni = 1, 712 = 2, the original integral is finite, but all terms from the hard region 
are infrared-singular and all terms from the soft region are ultraviolet-singular, even the 
leading-order contributions with j = and ji = j2 = 0. 

The contributions (2.8) and (2.9) can easily be evaluated. The hard-region integrals 
are straightforward by standard methods. They yield 

_ ,.2e .7E ( 2 •^^2-nl2-. ( ^(2 - ni - e) 



r(ni)r(n2) 



r(7ii2-2 + e + j)r(2-n2-e-j) 

j\T{A-n,,-2e-j) ^^'^^^ 
and can be summed up to the all-order hard contribution 

oo 

_p(/i) _ ^ jr(^) — ge7E g-j7rn.i2 (_p2 _ ^Q^2-ni2-e 
j=0 

r(ni2 - 2 + e) r(2 - ni - e) r(2 - ^2 - e) 



X 



r(ni)r(n2)r(4-ni2-2e) 



ni2-2 + e, ni2-3 + 2e; ^2 - 1 + e; j , (2.12) 

where 2F1 is the hypergeometric function and "— zO" indicates the side of the branch cut 
for the analytic continuation in the case p^ > 0. 

The soft-region contributions can e.g. be solved by tensor reduction (for j2 > 0). The 
integrals are only non-zero for even j2- We can identify j = ji + ^ as the number of 
additional powers of m?' /p^ compared to the leading order and rewrite 

00 00 

E^i. = E^^^- (2-13) 

ii,j2=o i=o 

An evaluation of the soft contributions (2.9) in closed form without the need for tensor 
reduction is shown in appendix A. 1.1. The result reads 

rfa+,)r(„.-2 + . 

j!r(2-.!i-j-i) ^ ' 
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and, summed up to all orders, 



pis) ^^p{ , ^ g 



2e pe7E ^,-i7™i2 ^^2-j2-n2-e 



- iO) 



j=0 



r(n2 -2 + 6) 

r(n2) 



m 



X 2-Fi ( ni, ni - 1 + e; 3 - n2 - e; ) . 



(2.15) 



The asymptotic expansion of the original integral (2.1) in powers of m /p is obtained 
by adding the contributions of the hard and the soft region. Which of the terms F^^^ and 

F^**^ are of the same order in /p^ depends on n2 (and e). For the special choice ni = 1, 
n2 = 2, we have 



Ah) 



1 f m 



is) 



\ p 

1 f m 



-p^ — iO 



e^7Er(l + e)r(l-e)r(-e) (2e)j 



r(l - 2e) 



p^ \ p 



/^2xe 



(eh 



(1 - e)j ■ 



(2.16) 



All terms -F^-^*^ are infrared-singular and all terms Pj^'^ are ultraviolet-singular, but only 



the leading-order terms Fq'^^ and Fq^' exhibit an explicit 1/e singularity: 

-p^ - iO 



As) 



pih) _ 1 

-p2 



1 , 

h In 

e 



+ 0{e) 



pih) 



2 fm^y 1 , 



F 



is) 



1 



1 



+ In 



F, 



is) 



1 fm'^V 1 



p^ \ p- 



+ 0(e), (2.17) 



where the column to the right is valid for j > 1. The 1/e poles cancel each other such that 
the complete expansion terms 



77. _ pih) , pis) 



(2.18) 



are all as finite in the limit e — )■ as the original integral: 



Fo = ^ In 



-p^ - zO 



1 /m^V 1 ^ 

-- — - + 0{e),j>l. 
pz ypZ J J 



(2.19) 



But remember that in Fq an infrared 1/e pole has been cancelled by an ultraviolet pole, 
which looks unnatural. We will come back to this point in the next section. 

The known exact result is reproduced by summing up all orders of the asymptotic 
expansion, either from (2.19) or by expanding the summed contributions (2.12) and (2.15) 
from the regions about e = 0: 



F = V F = — 

Z—^ J »j2 

i=o ^ 



In 



-p"^ — iO 



m 

+ ln( 1-^ 

pZ 



+ 0(e). 



(2.20) 
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For general ni and 722, the original integral can also be evaluated directly by standard 
methods, yielding 



p _ ,.2e .7E ^-^^n,, , 2.2-n,,-e r(2 - m - e) r(ni2 - 2 + e) 

r(n2)r(2-e) 



X2i^ifni,7ii2-2 + e;2-e;^^^^j . (2.21) 

By applying the transformation formula for the hyper geometric function which inverses 
its argument (see e.g. [14]), the sum of F^^^i (2.12) and F^**) (2.15) is exactly recovered 
from (2.21). Another check is obtained by evaluating the full integral with the help of 
a Mellin-Barnes representation and extracting the complete series of residues with rising 
powers of m? /p^. By doing so, the expansion terms (2.11) and (2.14) are reproduced. 

As a final remark for this section I want to emphasize that — although the evaluation 
of the integrals originating from the expansion by regions to all orders in the expansion 
parameter can be quite tedious — it is usually rather easy to extract just the leading- 
order contributions from all regions. This method is therefore particularly well suited for 
obtaining the leading order or the first few terms of an asymptotic expansion. 

2.2 Proof of the large-momentum expansion 

While the previous section has introduced the use of the expansion by region, we restart 
in this section from the original integral and transform it in a mathematically well-defined 
way until we arrive at the expanded expressions which have already been employed and 
evaluated in the previous section. 

We want to use the expansions of the hard (2.4) and soft (2.5) region which, obvi- 
ously, do not converge towards the original integrand throughout the complete integration 
domain. But we can divide the integration domain into a hard domain and a soft 
domain Ds, 

Dft = {/k G M'^ : |A;2| > A^} , = {A; G M"^ : |A;2| < A^} , (2.22) 

with some intermediate scale A^ chosen such that ^ A^ ^ These two domains 
are non-intersecting and cover the complete integration domain: 

Dhr\Ds = ^, DhLlD, = R'^. (2.23) 

When the loop momentum k is restricted to one of these domains, the corresponding 
expansion converges absolutely and the integrand is identical to its series expansion as long 
as the latter is summed to all orders: 



j 

I = T^-^)l = T^I for keDs. (2.24) 
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For the soft expansion, the summation index j symbohcally represents a proper combination 
of the indices ji and j2 in (2.5). The hard expansion r(^) only requires |/c^| ^> m? which 
is certainly fulfilled for k S D^- But the statement in (2.24) about the soft expansion 
is less trivial: While one of the conditions, <^ surely holds for k Ds, the other 
condition, |2A: - pi <^ can still be violated. We have to remember, though, that the 
expansions are always performed under the loop integral jDfc. By tensor reduction, each 
power of {k ■ p)^ in the numerator of the soft-region integrals (2.9) gives a contribution 
proportional to k'^p^, while odd powers of k ■ p vanish under the integration. This is still 
true if we restrict the integration by the Lorentz-invariant condition | A:^ | < . So we can 
safely count \2k - pi ~ \k'^p'^\^/'^ , and the condition \2k - pi ^ \p'^\ holds under the integral 
over the soft domain Dg. 

Alternatively, for ensuring the convergence of within D,, one may choose a refer- 
ence frame in which either the zero-component or the spatial components of the vector p 
vanish, depending on the sign of p^, and define the boundaries of Dh and Ds appropri- 
ately in this reference frame. Or, for p'^ < 0, one may perform a Wick rotation and define 
the boundaries of the domains as relations between positive-definite norms of Euclidean 
vectors. 

A consequence of the absolute convergence of the expansions is that the summation 
commutes with the integration and can safely be pulled out of the integral if the latter is 
restricted to the corresponding domain: 

j ^kI = Y^ j BkT^'^h , j BkI = Y^ j BkT^'h. (2.25) 

After these preliminaries we can start transforming the original integral by splitting the 
integration into the two domains and performing the appropriate expansions in each of 
them: 

F = jBkI= j Bkl+ j BkI = Y^ j Bkrj'^h + Y^ j Bkrfl. (2.26) 

fceDh. keDs ^ keDh ^ k&Ds 

With dimensional regularization at hand, we can also perform the integration of each 
expanded term over the complete integration domain, but we have to compensate for this 
by subtracting the surplus from the added domain: 

BkT^'^h = [BkT^^h- I BkT^^h. 



3 

k&Dh k&D. 



BkTf'I = J BkTf'I- J BkTf'I. (2.27) 

keDs keDh 

Without the indication of a restriction, the integrations are understood as being performed 
over the complete integration domain W^. Relation (2.25) also holds if the integrand I is 
not the original one, but a term from a previous expansion, i.e. each expansion can be 
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applied to any integrand if the integral is restricted to the corresponding domain: 

j BkT^'h = Y^ j mr^^hfi, j BkTj^h = Y^ j mT'>'h'-^h. (2.28) 

Let us have a look at the newly generated double expansions. The order in which the hard 
and soft expansions are applied is irrelevant because the doubly expanded integrand is the 
same in both cases (written here with two separate indices for the soft expansion to be 
specific) : 

(/i) (s) ^ (s) {h) J ^ {n2)i (ni)j,, (m^)' (-1)^1^ {2k ■ 

i iiJ2 iiJ2 i j^l j^ljg! (p2)ni+ji2 (^2)n2+i-ii ' \ ■ ^) 

In such cases of commuting expansions, we label multiple expansions by a comma-separated 
list in the round brackets: 

J,{h)rp{s) _ rp{s)rp{h) _ rp{h,s) 

After an appropriate relabelling of the summation indices, the two contributions with 
double expansions can be added together, 

EE / D^^S"^^+EE / D^7;(;-)/=j;|DA.rif)/, (2.31) 

arriving at an integral over the complete integration domain. The non-trivial point here 
is that we have to exchange the order of the two summations in one of the contributions. 
While e.g. in the first term with k restricted to the soft domain Ds, the summation 
of the soft expansion is absolutely convergent, we cannot easily claim convergence for 
the summation of the hard expansion when k £ Dg. However, we are not summing 
expanded integrands here, but integrals. And the only scale involved in the integrals 
over the doubly-expanded integrand T^^-'^^ I (2.29) originates from the boundary of the 
integration domain, as all occurrences of the momentum p in the scalar products k ■ p m. 
the numerator can be pulled out of the integral. In fact, for dimensional reasons, we know 
that 



because the only dimensionful parameter in the definition of the domain Ds is A^. By 
the same reasoning, (2.32) holds if the integral is restricted to A; G instead. As the 
boundary has been chosen to obey m? <C A^ <C both the summations over i and over 
ji, j2 in (2.32) converge absolutely, and their order can be exchanged at will. 
We are now able to collect all pieces contributing to the integral F. Writing 

^ = jm rfh , = jm rfi , F^y^ = fm T^yh (2.33) 
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for the individual contributions and 



p{h) ^ p{h) ^ p{s) ^ ^ p{s) ^ pih,s) ^ p(h,s) (2.34) 

3 3 hj 

for the summed- up series, we obtain 

F = f(^) + f(^) - F'-''''^ (2.35) 

for the original integral after the above transformations. Note that all integrations involved 
in (2.35) are performed over the whole integration domain MJ^. So all restrictions to the 
two individual domains and Dg drop out and the final terms in (2.35) are individually 
independent of the separating scale A^. Thus the exact position of the boundary between 
the domains is irrelevant, and we could have defined the domains e.g. in the following way: 

Dh = {k £ R'^ : \k'^\ ^ m'^] , Ds = {k £ R'^ : \k'^\ < m'^] , (2.36) 

where "<" is understood as the negation of "S>" , such that r\Ds = % and VJDg = R'^ 
hold. In later examples we will not introduce specific boundaries between convergence 
domains, but use rather "sloppy" specifications of the domains as above. It is understood, 
however, that exact positions of the boundaries exist and could be specified if needed. 

The first two terms in the final identity (2.35) correspond exactly to the contributions 
from the hard (2.8) and soft (2.9) regions prescribed by the expansion by regions and 
evaluated in the previous section. But now we have obtained a third term, subtracted from 
the first two. This additional overlap contribution F^^'^^ is absent in the recipe formulated 
in section 1. Let us have a look at its terms: 

pih,s) _ [r.j,r^h,s) _ in,), (m),,, jm^y f {2k ■ _ , . 

These are scaleless integrals, which must consistently be set to zero when using dimen- 
sional regularization. In fact, each of the integrals in (2.37) can be transformed by tensor 
reduction into an integral /DA; (/c^)~" with some power n, and these massless tadpole inte- 
grals exhibit both ultraviolet and infrared singularities in such a way that they cancel each 
other. (Individual parts of the integration have different convergence domains in the com- 
plex plane of the space-time dimension d, but analytic continuation permits to combine 
the pieces.) The integral 

DA; 1 1 , , 

712X2 = = (2-38) 

[k-^y euv eiR 

is the only case in this class of integrals where ultraviolet poles 1/euv or infrared poles 
1/eiR appear. In the special case ni = 1, n2 = 2, where the original integral F is finite, the 
integrals of the overlap contribution (2.37) exhibit exactly those ultraviolet singularities to 
cancel the ones in the soft contribution F^^^ and those infrared singularities to cancel the 
ones in the hard contribution F^^\ Although the overlap contribution F^^''^^ is scaleless and 
vanishes, it is this term which makes the complete result (2.35) separately ultraviolet-finite 
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and infrared-finite. One can check explicitly (see appendix A. 1.2) that the summed-up 
overlap contributions (2.37) with ultraviolet and infrared 1/e poles separated cancel the 
corresponding poles in (2.17). 

Note that the overlap contribution terms (2.37) are the same as the doubly expanded 
terms arising in an expansion by subgraphs [4]. 

Having checked that the overlap contribution F^^'^^ is scaleless and vanishes, the orig- 
inal integral is reproduced by 

F = F^^) + F^'^ . (2.39) 

This is exactly the sum of contributions which has been evaluated in section 2.1, where 
we have assumed that we need these two regions (hard and soft) and evaluated them 
according to the recipe of the expansion by regions. Now we have obtained the same 
answer by mathematically transforming the original integral. Let us recapitulate what we 
had to check on our way: 

• For the two regions we had to find domains {D^ and Dg) where their expansions 
converge absolutely. These domains have to be non-intersecting (0^ H Z?<j = 0) and 
cover the complete integration domain (Dh. U Ds = M."^). 

• In the double expansion the order of the two expansions has to be irrelevant ( "com- 
muting expansions"). 

• The overlap contribution from the double expansion involves only scaleless integrals. 

These three points had to be proven explicitly for the example integral. The rest of the 
transformations used in this section is general and applies to any other integral with a set 
of regions and domains obeying analogous conditions. We will work this out in the next 
section. 

Note that we did not have to evaluate any of the integrals in F^^\ F^^^ or in 
order to prove the identity (2.35). It is sufficient to study the expansions at the integrand 
level. And even for the form (2.39), where the scaleless overlap contribution has been 
dropped, a look at the expanded integrand has been enough (although in other cases it can 
be more involved to show that the overlap contributions are scaleless). 

Remember how important it is within the framework of the expansion by regions that 
scaleless integrals can be set to zero. In our example it is dimensional regularization which 
regularizes scaleless integrals in such a way that they vanish. In some cases (see in particular 
the examples in sections 6 and 7) this is not sufficient, and we have to use analytic regu- 
larization as well. In the absence of such nicely behaved regularization schemes, however, 
interesting patterns appear (see section 7.3 and appendix B) where overlap contributions 
play an important role. 

3 General formalism with commuting expansions 

In this section the proof of section 2.2 is generalized to arbitrary integrals, with some 
restrictions. Consider the following situation: 
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We want to expand the integral 

'Bkl, (3.1) 



where the integrand / is integrated over the integration domain D. This can be a 
one- loop integral (with D = M'^), a multi-loop integral {D = M" '^) or any arbitrary 
integral. 

• We have identified a set 

R = {xi,. . . ,xn} (3.2) 
of N regions Xj. Each region x is characterized by an expansion 

rpix) = Y^T^-) (3.3) 

j 

which, when applied to the integrand, replaces the latter by a series of expanded 
terms. The summation index j can also represent a set of indices ji, j2, . . ., but we 
only write one index per expansion. These expansion operators also have to be defined 
when they are applied to terms resulting from previous expansions (when multiple 
expansions are generated). In such cases it may happen that a certain expansion 
is an identity transformation and the set of summation indices represented by j is 
empty. 

• For each region x there is a domain C D such that the expansion T^^^ converges 
absolutely when the integration variable is restricted to k £ D^- 

Let us assume that the regions, expansions and domains fulfill the following conditions: 

1. The domains are non-intersecting and cover the complete integration domain: 

D^r\D^, ={byx^x' , [jD^ = D. (3.4) 

x&R 

2. All expansions commute with each other: 

r(^)r(^') = r(^')r(^') Vx, x' g i? . (3.5) 

3. The original integral itself and all integrals over expanded terms — whether restricted 
to some convergence domain or not — are regularized. 

4. The series expansions converge absolutely (or are properly regularized) even 
when the expanded terms are integrated over the whole integration domain D instead 
of just their convergence domain Dx- 

Condition 1 ensures that any integral over the complete integration domain can be split 
into integrals restricted to the domains Dx'- 

m= j m = Y^ j ^k. (3.6) 

k&D k&D:, 
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The consequence of this condition is that we might have to invent "auxihary" regions 
(which do not contribute to the final result) in order to cover the complete integration 
domain D with the convergence domains D^. 

Condition 2 is to be understood at the operator level of the expansions: Whatever 
integrand the double expansion T^^^T^^ ^ is applied to, whether the original integrand / or 
a term resulting from previous expansions, the order in which the expansions are performed 
must be irrelevant, i.e. in all cases the same (multiple) series of doubly expanded integrand 
terms is established. This condition, however, cannot be fulfilled in all cases where the 
expansion by regions has been applied successfully. In section 5 this restriction will be 
relaxed and the treatment of non-commuting expansions will be presented. 

Condition 3 implies that we use a regularization prescription which provides a math- 
ematically well-defined meaning to all integrals occurring in the formalism described in 
this section. Usually this is dimensional regularization, eventually combined with analytic 
regularization, but other schemes are possible. 

Finally condition 4 requires a certain mechanism that makes series expansions converge 
even outside their convergence domains. For loop integrals with dimensional regularization 
this is usually the case. More generally, this mostly works for integrals where the boundaries 
either lie at zero or at infinity such that they do not introduce new scales. Then the 
integration over the complete domain Z), although formally divergent, is regularized and 
determined only by the scaleful parameters within the convergence domain, keeping the 
series expansions as convergent as with integrals restricted to the convergence domains Dx- 
It is possible to apply the expansion by regions to other integrals, e.g. with domains D 
involving finite boundaries, where condition 4 is possibly violated. But then the summation 
of the series expansions has to be done with care, and it might be necessary to combine 
certain terms which are individually divergent. This is a subtle issue in the expansion by 
region and in the formalism presented here. Its consequences will be pointed out at the 
relevant steps later in this section, and appendix B presents an example illustrating this 
behaviour. 

Let us introduce the notations to be used in this section. They have partially already 
been defined for the example in section 2. Multiple expansions (replacing some integrand 
first by its expansion terms according to one region and repeating this step with the re- 
sulting terms for other regions) are denoted by 

rj,{xi,X2,...) _ rp{xi)rp{X2) _ _ _ rp{xi ,X2 ,. . .) — rp{xi,X2,...) 

if the expansions are commuting, i.e. their order is irrelevant. The j-th order expanded 
integral according to the region x is denoted by 

Fj"^^ = J BhT^'^h, (3.8) 

and its summation to all orders by 

3 3 
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where the integrals are performed over the complete integration domain D. Analogous 
notations are used for multiple expansions: 



p{xi,X2,...) _ j^j^rp(xi,X2,...) J p{xi,X2,...) ^ pix:i,X2,...) (3 10) 

The restriction of an integration to a domain is indicated by a lower index in square 
brackets: 

F'[x]= j Dfc^, 4"'n= / DA:Tj";-)/, etc. (3.11) 

k&Dx keDx 

If the integration is performed over the combination of several domains, we write 

n n 

F[xi + ... + X„] = '^^[X^] 1 ^j,..'lxi + ...+X„] = '^^j,..'.lx,] ■ (3.12) 

1 = 1 1 = 1 

The absolute convergence of each expansion T^^^ within the corresponding domain Dx 
implies that we can safely replace any integrand /' (the original integrand / or the result 
of previous expansions) by its expanded series if the integration variable is restricted to 
k £ Dx. 

I' = T^^)l' = ^Tfl' fovkeDx. (3.13) 

j 

Absolute convergence also implies that we can pull the summation of such an expansion 
out of an adequately restricted integral, using the notation of (3.8)-(3.11): 

j j 

Now we can start with the original integral (3.1) and split the integration into the 
N domains corresponding to the regions, according to (3.6). In each domain we replace 
the integral by its series expansion according to (3.14): 

xGR a;Gi? 

The right-hand side of (3.15) involves a sum of N series expansions with each integral 
restricted to the corresponding convergence domain. This is a special case (with n = 1) of 
the expression 

E ' (3-16) 

{x{,...,x'„}cR 

where the sum runs over all subsets of n distinct regions out of the regions in i? (1 < n < 
N). Each integrand is multiply expanded according to these n regions, and the integrals 
are performed over the combination of the n corresponding domains. Let us postpone for 
a few lines the question whether the expression (3.16) is a convergent series expansion, 
despite the fact that (for n > 1) the integrations are performed over larger domains than 
the convergence domain of each of the n individual expansions. 
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If n < A^, i.e. if the integrations in (3.16) are not performed over the complete inte- 
gration domain D yet, the regularization of the integrals (see condition 3 above) allows 
to extend all these integrations to k £ D when compensating for this by subtracting the 
integrations over the additional domains: 



[x\+...+x'„] ~ jl,...,j„ [x'^- 



Jl.-Jn x'^^^(^R\{x'^,...,x'^} 



The subtraction terms are integrals performed over one domain D^t ^ each (where x'^^^ 
runs over all regions which are absent in the subset {x^, . . . , x^}). These subtraction terms 
can be replaced by their expansions according to (3.14): 

jn + l 

where we have already used condition 2 above that the expansions commute. Now we sum 
the individual terms in (3.17) separately and write 

<+ie-R\K,...,<} 

This is a non-trivial step. Even if the complete expression (3.17) is a finite series expansion, 
this does not necessarily mean that all summations in (3.19) are individually convergent. 
Depending on the parameters involved, especially the boundaries of the complete domain D 
and the boundaries between the individual domains D^, we might have to regularize the 
summations in (3.17). We may e.g. think of truncating the summations by some upper 
limit (ji < jmaxVi), thus limiting the accuracy of the expanded expressions, but dealing 
only with finite sums. This truncation is removed (by jmax oo) only in the end when 
the summations have been combined into convergent ones (cf. condition 4 above). 

For reproducing (3.16) we finally we have to sum (3.19) over all subsets of n regions 
{x'l, . . . ,x'^}. The subtraction terms yield summations over subsets of {n + 1) regions, 
where each term appears {n + 1) times with different integration domains D^' ^^: 

n+l 

{x[,...,x'^}cR x'^^,eR\{x[,...,x'„} {x[,...,x'^^,}cR i=l 

This requires not only that the expansion commute (condition 2), but also that the order 
of the series summations in (3.20) can be exchanged, which in turn requires their absolute 
convergence. For the example presented in section 2.2 this has been shown explicitly. 
I cannot provide a rigorous proof for this convergence issue in the general case treated 
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here, but I am convinced that even divergent series in intermediate steps of this derivation 
(which then have to be regularized) are not problematic when condition 4 ensures the 
convergence of those sums which remain in the final result. Note that this problem only 
arises when the series expansions are summed up to all orders. If an approximation with 
limited accuracy is sought and the series expansions are truncated at some finite order of 
the expansion parameter, then convergence problems of individual (intermediary or final) 
terms are absent. 

Combining all terms from (3.19) and (3.20), the expression (3.16) yields 

K,...,<}ci? {4,...,<}ci? K.-,<+i}ci? 

The first term on the right-hand side consists of integrals performed over the complete 
domain D. We want to keep such terms for the final result. The second term is exactly the 
same as the one on the left-hand side, but with n replaced by (n + 1) and with opposite 
sign. Thus (3.21) represents a recursion formula which can be iterated from n = 1, as 
in (3.15), up to n = — 1. This allows us to write the original integral in the following 
form: 

_p = ^ _ y^ .4) _)_..._ y^ 

x£R {x{,x'2}cR {x[,...,x'„}cR 

+ ...-(-1)^f(^1'-'^^). (3.22) 

This is the master identity for the expansion by regions in the formalism presented in this 
section. It involves only integrations over the complete domain D, so at least the series 
expansions in this final result are all convergent if condition 4 holds. 

The master identity (3.22) involves single and multiple expansions, according to the 
N regions xi, . . . ,xjsf and all of their combinations, with alternating signs. The recipe 
for the expansion by regions presented in section 1 only knows about the first term on 
the right-hand side of (3.22), where a single-expanded integrand according to each region 
is integrated over the complete domain. This means that all other terms with multiple 
expansions must vanish in "normal" situations for the recipe to be valid. 

Indeed, if the regularization of a loop integral and the regions are chosen properly, 
then each of the terms Fj^^ (i.e. the single-expanded terms present in the known recipe) 
is a single-scale integral yielding a homogeneous function of the expansion parameter with 
a unique scaling^ (cf. section 2.1). Every further expansion of such a single-scale integral 
according to a different region (which would yield a different scaling with the expansion 
parameter) makes the integral scaleless such that it is set to zero. This is why usually 
the terms with multiple expansions do not contribute to the asymptotic expansion of loop 
integrals. 

^This means that each term in the series of expanded integrals depends on the expansion parameter by a 
simple power which is different for each region. If two regions share the same dependence on the expansion 
parameter, then their overlap contribution does not have to be scaleless. Note that single-scale integrals as 
understood above may exhibit a non-trivial dependence on additional parameters or ratios of 0(1) which 
are considered neither small nor large in the expansion. 
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But the identity (3.22) is more general: It is independent of the chosen regularization 
scheme, as long as all individual terms are mathematically well-defined (as required by 
conditions 3 and 4). The identity (3.22) can also be applied to other types of integrals or to 
loop integrals which use other regularizations than the standard dimensional and analytic 
ones. In such cases the overlap contributions, i.e. the terms in (3.22) with expansions 
according to more than one region, may become relevant. 

While these overlap contributions arise naturally out of the formalism presented here, 
they have already been noted in the context of effective theories and are called "zero-bin 
subtractions" there (see e.g. [11, 12]). The overlap contributions in (3.22) have exactly the 
same form of multiple expansions as what is — usually only to leading order — introduced 
in the literature under the name zero-bin subtractions. The identity (3.22) clarifies the 
whole picture of subtractions which are needed in the general case. 

One example where overlap contributions (or zero-bin subtractions) are relevant is 
provided in [12] where A-regulators are introduced which push all propagator denominators 
artificially off-shell by some amount. These A-regulators introduce new scales into the 
integrals. Therefore the contributions from each region are not homogeneous functions 
of the expansion parameter and the overlap contributions are not scaleless. Note that the 
authors of [12] only consider the overlap of each collinear region n, n with the soft region s, 
in the language of this section and The identity (3.22) would have told them 

that the full result is (assuming that this set of regions is complete) 

F = F^""^ + F^""^ + F^'^^ — F^"^'^^ — F^"''^^ — F*^"'"-* -|- F*'"'"'"*-' . (3.23) 

But the soft expansion of their integral is identical to the combination of the two collinear 
expansions, such that F'^'^ = F^"'") = F^"'^) = F^^''^ = F^"'"'^). Cancelhng the last two 
terms in (3.23), the result can be written as 

F = ^F*^"^ - F*^"'"*)^ -I- ^F^") - F^"'"*^^ -I- F*^**) , (3.24) 

as they did, or by omitting the (irrelevant) soft region: 

F = f(") + f(") - f("'^) . (3.25) 

The master identity (3.22) is an exact relation for the original integral F. It involves 
the summation of the single and multiple series expansions to all orders according to the 
definitions (3.9) and (3.10). If instead only the leading terms of the right-hand side in (3.22) 
are taken into account, then an approximation for the integral F on the left-hand side is 
obtained. The expansion parameter by whose powers higher-order terms are suppressed is 
related to the parameter hierarchies exploited in the expansions T^^\ While intermediate 
expressions in the derivation of this formula involve the boundaries between the convergence 
domains, the final result (3.22) is independent of these boundaries, and so are the series 
expansions in the master identity. 

A leading-order asymptotic expansion of the integral F can also be obtained directly 
without ever touching infinite series expansions. If, in the steps (3.15) and (3.18) above. 
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the integrands are simply replaced by their leading-order expansion terms according to the 
regions x and x'^_^_^, respectively, 

^^Z^^oN' V..,oK+i] ^^o,...,o,o[<+i] ' ^'^■^'^) 

x&R 

then higher-order terms are neglected which are suppressed either by powers of the expan- 
sion parameter or by other small parameter ratios involving the boundaries between the 
convergence domains. Remember that these expansions are only introduced when they are 
absolutely convergent because the integration variable is restricted to the corresponding 
domain, or D^./ , respectively. In the course of the derivation above, with all sum- 
mations replaced by their leading terms, the contributions are combined such that any 
dependence on the boundaries cancels out. Finally the leading-order approximation 



x&R {x'^,x'^}(lR {x'-^,...,x'„}(lR 

+ ..._ (3.27) 



is obtained, which reproduces the original integral F up to terms suppressed by powers of 
the expansion parameter. The leading contributions from the different regions may start 
at different powers of the expansion parameter such that only some of the terms in (3.27) 
are actually leading-order contributions while others are suppressed. 

The leading-order expression (3.27), although it may be derived from the all-order 
result (3.22), has a value of its own because it still holds when the validity of condition 3 
above cannot be verified for higher-order terms or when condition 4 is violated because 
the summations do not converge individually. Appendix B illustrates such a behaviour 
of the expansion by regions with an example involving a finite integration boundary and 
non-converging series expansions. 



4 Example: threshold expansion 

Before generalizing the formalism of section 3 to non-commuting expansions, let us apply it 
to another example. The threshold expansion of the one-loop three-point integral presented 
here is also the first example treated in the original paper [1]. It is illustrated in figure 2. 
We choose the loop-momentum parametrization which the authors of [1] only use for the 
soft/ultrasoft region because this parametrization is better adapted for what we want to 
demonstrate here: 

F = /d/e / , with / = hhh and 

1 1 



((|+p + A:)2-m2)"^ {B +q-k + 2p-kY^' 



^2 ((|_p_^)2_^2)-2 {k^-q.k + 2p-kY^' (fc2)-3- (^-^^ 
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Figure 2. Loop integral for the threshold expansion. 



We will only evaluate this integral for ni = n2 = n-s = 1, but we need the general propagator 
powers for the analytic regularization of some contributions. In the expressions (4.1) for 
the propagators the on-shell conditions (| ip)^ = m'^ have been used. These also imply 
q-p = and = — We are interested in the threshold regime q'^ ~ (2m)^, where 

q"^ > l/l . (4.2) 

Let us choose as a specific reference frame the centre-of-mass system of the momentum q, 
where 

{qn = {Qo,0) and (p^) = (0,p). (4.3) 

The propagators then read 

A ' ■ ' 



{kl - /c2 + q^kQ - 2p- A:)"i (/eg - k"^ - qoko - 2p- /c)"2 

^3 = ^ , (4.4) 

and we are looking for an expansion in powers of \p\/qQ <^ 1. The loop-momentum com- 
ponents /cq and k are multiplied with prefactors of different orders of magnitude in the 
propagators Ii and l2- Thus it is natural that we also get a region which distinguishes 
between /cq and k. The three regions needed for this example are 

• the hard region {h), characterized by /cq ~ A; ~ q^, with the expansion 

rW/i2 = y^ , rW/3 = /3, (4.5) 

U ^- {ki-k^±q,k,r^^-^r ^ ^ 

converging absolutely within = {k D : \k\ ^ \p\ V \ko\ 3> \p\}, 

• the soft region (s), characterized by /cq ~ A; ~ p, with the expansion 

-'''-'1,2- > YT , yri, + ' -t -'3--'3, (,4.bj 

Jl J2J3=0 

converging absolutely within Dg = {k £ D : |A;|<|A;o|<|p|}, 
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and the potential region (p), characterized by ko ~ p'^/qo and k ^ p, with the 
expansion 



j=0 



f- {-P ± qoko - 2p- fe)"i.2+i 



Vi!^ JzM)^ (47) 

converging absolutely within Dp = {^k ^ D : |^o|^|^|^|p|}) 

where = M*^ is the complete integration domain. We do not have to specify the exact 
positions of the boundaries between the convergence domains, and the relation "<" is 
understood as the negation of like in (2.36). So condition 1 of the formalism in 

section 3 holds: 

Dhr^Ds = Dhr^Dp = Dsr^Dp = ^, DhUDsUDp = D. (4.8) 

All expansions commute with each other (condition 2), and the multiple expansions read: 

-P ± goA;o)"'-2+^'''' ' 

rp{h,S,p)j,^ = r(P)/3 . (4.9) 

Appendix A. 2.1 demonstrates that the expansions (4.5)-(4.7) converge absolutely when 
the integration variable k is restricted to the corresponding domain. The domains D^, 
Dg and Dp have been chosen as large as possible within the convergence domains of the 
expansions. They contain parts which do not correspond to the scaling of the loop momen- 
tum components specified for each region. We may e.g. have \k\ <C |A:o| within Dg, which 
contradicts k ^ k^ ^ p. But by choosing such enlarged domains we avoid the introduction 
of additional, artificial regions for covering the complete integration domain. 
The expanded integrals read 



{h,s) : 




= tW/i,2, 


{Kp) : 




31,32=0 


(s,p) : 






{h,s,p) : 







Bk {2p- ky^^ 



(A;2 -P + qok^Y^+h {kl -P- qokoY^+n (yfcg - Py^ ' 

p{s) ^ (^l)jl23 (^2)j456 

Jl,...J6=0 



Dfc {-kiy^* (P)J25 (2p • a!)J36 



{qoko + iO)"i+Ji23 (-go^o + i0)"2+i456 (kl - k'^Y^ 
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pip) ^ (ni)ji {n2)j2 (n3)j3 

rf n il!j2!j3! 



(-/c2 + q^ko - 2p- kY^+n (-k^ - q^ko - 2p- /c)'"2+i2 (_/j2)n3+i3 



p{h,s) _ p{s) ^ 

p{h,p) ^ (ra2)j34 (n3)jg 



J1,...J5=0 



(ni)j,23 ('^2)i456 ('^3)i7 



n Jl!---J7! 

J1,...J7=0 



(gO^O + iO)'"l+Jl2^ (-gO^O + i0)"2+i456 (_^2)n3+j7 

p{h,s,p) _ p{s,p) ^ (4-10) 

As we will see in a moment when evaluating the expressions in (4.10), the integrals there 
are well-defined through dimensional and analytic regularization, and the summations are 
absolutely convergent. So also conditions 3 and 4 of section 3 hold, and the original 
integral (4.1) can be expressed through the master identity (3.22): 

p — p{h) _|_ p{s) _|_ pip) _ pih,s) _ p{h,p) _ pis,p) _|_ p{h,s,p) _ (4-11) 

Now, independent of what the contributions F^*^ and F^^'^^ are, we can see that they 
drop out because the hard expansion T^^^ does not change the integrand if the latter has 
been expanded via T^*) before, i.e. F^'^'^^ = F^^^ and F^'*'*'^) = F^^'^), and therefore 

(^(■5) _ F^'^'^^) - (F*^*'^^ - f'^'*''*'^^) = 0. (4-12) 

So all terms including the soft expansion in (4.11) do not contribute to the result. 

Examining these terms nevertheless, we recognize them as scaleless contributions: Scal- 
ing the loop momentum by kQ — )■ A/cq and k ^ Xkin the integrals of F^^^ and F'-^'P) in (4.10), 
we notice that each of these integrals is identical to itself times A^~"'i2-2n.3+ii245-2e^ Within 
dimensional regularization (e Z) and analytic regularization (n^ Z), this factor is dif- 
ferent from 1, so the integrals are scaleless, they either vanish or diverge. Without analytic 
regularization, the integrals of F^**) and F^'^'P^ are regularized by the (3 — 2e)-dimensional k- 
integration, but the one-dimensional /cQ-integration can still be divergent when it is consid- 
ered separately and evaluated before the fe-integration. In particular, the /co-integration is 
then singular because the integration contour is pinched between two poles both at A;o = 0, 
but on different sides of the contour. ^ So, strictly speaking, we need analytic regularization 



^In [1] the authors argue that the pinching of poles must be ignored in the soft region because these 
poles have already been taken into account through the potential region. In the formalism developed here 
we cannot use this argument because we have to take the integral as it arises from the expansions, and all 
possible double-counting is eliminated via the subtractions of the overlap contributions. 
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(through non-integer powers of the first two propagators) or some other additional regu- 
larization here in order to make the integrals well-defined. Then the contributions F^^\ 
^ jp{s,p) g^j^jj p{h,s,p) g^j^g scaleless and must be set to zero. 

A similar argument shows that the contribution F^^''^^ in (4.10) is scaleless as well: 
Scaling the loop momentum components by kQ — )■ X^ko and k — >■ Xk, each of these integrals 
is found to be identical to itself times A^-2"i23+2ji35-i24-2e^ ^i^^ p(h,p) vanishes. 

Finally we obtain 

The evaluation of these two remaining contributions is sketched in the appendices A. 2. 2 
and A. 2. 3. For n\ = n2 = nj, = 1 the results read [1] 



q2 J ji (1 + 2e + 2j) ' 

Note that the potential region only has a leading-order contribution, its higher-order in- 
tegrals vanish, cf. appendix A. 2. 3. The series expansion of the hard contribution F^^^ is 
absolutely convergent for \p'^\ <C as required by condition 4 in section 3. The hard 
contribution can be summed up and expressed through a hyper geometric function: 

The original integral (4.1) can alternatively be solved by introducing a Mellin-Barnes 
representation. For ni = 77-2 = ^^3 = 1 the Mellin-Barnes integral reads 



q^\q^J e J 2nT \ q^ ) -I - e - z 

—ioo 

where the pole at z = —■^ — e lies to the right of the integration contour. An asymptotic 
expansion in powers of jq^ is obtained by closing the 2;-contour to the right. The residues 
of the poles of r(— z) reproduce the hard-region expansion F'-'^^ in (4.14), whereas the 
residue at 2; = — ^ — e yields the potential contribution F^^^ . 

The Mellin-Barnes representation (4.16) also permits to perform an asymptotic ex- 
pansion in the opposite case, \]?\ S> (7^, by closing the contour to the left. Here only one 
series of residues, from the poles of r(l + 6 + 2:), contributes and is expressed through a 
single hypergeometric function: 

Through analytic continuation of this result to ^ \ by inversing the argument of the 
hypergeometric function (see e.g. [14]), the sum of (4.15) and (4.14) is recovered. 
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5 Formalism for non-commuting expansions 



Let us return to the general case. When introducing the general formalism in section 3, 
we required (condition 2) that all expansions commute with each other. As we will see in 
examples, this condition cannot always be fulfilled by a proper choice of the regions. In 
the following paragraphs a generalized formalism is developed which relaxes this condition 
to some extent. 

We start with the same situation as described at the beginning of section 3: The 
integral F = jDk I with integration domain D shall be expanded. We have a set R of 
N regions, R = {xi, . . . ,xj\[}. Each region x is characterized by an expansion 

Tj which converges absolutely within the domain D^. Conditions 1, 3 and 4 hold, i.e. 
the domains are non- intersecting and cover the complete integration domain D, all integrals 
over expanded terms are regularized, and the series expansions F^^'^ (with integrals over 
the complete domain D) converge absolutely. 

Condition 2 is relaxed to the new condition 2a: 

2a. All expansions corresponding to regions within some subset Rc C R commute with 
each other and with expansions of any other region in R: 

rp{x)rp(x') ^ ri.(x')j.{x) G , x' ^ R . (5.1) 

In other words, two expansions can be interchanged if at least one of them belongs to a 
region from the subset Rc. Without loss of generality let the subset Rc contain the first 
A^c regions from R: 

Rc = {xi,...,xnJ C R, 0<Nc<N. (5.2) 

Let 

Rnc = R\Rc = {xn,+i,---,xn} (5.3) 

be the subset of regions with non-commuting expansions. Then two expansions are non- 
commuting only if their corresponding regions both belong to R^c- We do not specify how 
small or large the set R^c can be within R, but the formalism developed in the following 
will only provide useful statements if there are still regions left within Rc. Obviously the 
case Ac = A — 1 is equivalent to Ac = A because a single region within R^c would still 
provide an expansion commuting with all others, as specified in condition 2a above. 

The notations introduced in section 3 are used here as well, with one addition: Ac- 
cording to (3.7), a multiple expansion x^^^'^'^''"'^ = V ■ ■ rpixux2,...} jj^pj^gg ^j^^^ these 

(x) 

expansions T^, commute with each other, i.e. at most one of them is a non-commuting 
expansion of a region from R^c- Whenever two non-commuting expansions are applied 
successively, the order of their application has to be specified. We define 



„(a:2<— x'j) _ ^{x'2) rp{x'i) rp{x'2<—x'-,) _ \ ^ rp{x'2<—x[) I > TJ / ri 4\ 

h d2 
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as the operator which first expands according to the region x'l, then according to the 
region indicated by the arrow in the superscript. Such a pair of non-commuting expan- 
sions may be combined with further commuting ones corresponding to regions from Rc. 
They are specified in a comma-separated fist because the order of their apphcation with 
respect to each other and to the pair of non-commuting expansions is irrelevant: 

„{a;2<— XjjXg,...) — (X2) rp{x'i) rp{x':i) rpix'^) rpix2) rp{x'i) I /- D ( r. r,\ 

^i2,il J3,... - ^32 ^jl ^js "' - -^h "' ^32 ' X3, . . . /tc • l^-fJj 

Exactly as in section 3, we start by splitting the integration into the domains and 
expanding each restricted integral accordingly, see (3.15): 

x&R 

The right-hand side of (5.6) is a special case (for n = 1) of the expression 

(iJc + l> 



{x'^,...,x'J<zR 

where each integrand is multiply expanded according to n regions and integrated over the 
combination of the n corresponding domains. The sum runs over subsets of n distinct 
regions, as in section 3. But the superscript {Rc + 1) indicates that the sum is restricted to 
such subsets which contain at most one region from i?nc with a non-commuting expansion: 

(-Rc+l) 

E - E - E ■ (5-8) 

{x;,...,<}C-R {x'^,...,x'^}cR, {a:i,...,<_i}cRc, 

{x{,...,x'Jr\R^,=% or {x} x;efl\K,...,<_ J 

This restriction ensures that all expansions T^^i\ . . . , T^^"^ in (5.7) commute with each 
other and that the multiple expansion can be written in the form (5.7). Obviously such 
a restricted sum over subsets of regions excludes any pair of regions from /?nc with non- 
commuting expansions. 

Following the steps of section 3 in (3.17)-(3.19), we extend the integrations in (5.7) to 
the complete integration domain, compensating for this by subtraction terms: 

Note that, in contrast to (3.19), the subtraction terms in (5.9) have not yet been expanded 
according to the additional region x'^^^. We can perform this expansion T^^"+i\ but we 
have to distinguish two cases here: If {x[, . . . , C Rc or x^^]^ £ Rc, then all expansions 
T^^'i), . . . ,T^^"+i^ commute with each other and we proceed as in section 3: 

K,.. K,.. x;,x;,+,)_ (5.10) 
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But if 

x'n+i S -^nc and also {x'l, . . . , x^} already contains one region from i?nc (l^t us assume 
without loss of generality that x[ G Rnc and {xg, . . . , x'^} C Rc), then 

because the expansion T^^'i^ has been performed before T^^"+^\ Summing these terms 
according to (5.7) and (5.9), we get two contributions, one in analogy to (3.20) in section 3, 
the other involving pairs of non-commuting expansions: 

{Rc+l) 

EST pK,...,<) 

{x[,...,x'„}cR x'^_^-^eR\{x[,...,x'J 
{Rc+l) n+l 

{a:'l,-,<+i}Ci?^ i=l ^ x[eRnc x^giinc, {xJj,...,<+i}Ciic 

77i(^'l'---'^'n+l) 
[xi + ...+<+i] 



(5.12) 



We arrive at the following recursion formula for the expression (5.7): 

{Rc+l) {Rc+l) {Rc+l) 

Ax'„...,x'J _ r^{x[,...,x'„) _ V- ^W.-.<+i) 

^[x[ + ...+x'J - 2^ ^ ^[x[+...+x'„^,] 

{x[,...,x'„}cR {x{,...,x'jGR {xi,...,x'}cR 



Z — ^ [X2\ 



The first two terms on the right-hand side of (5.13) look similar to the ones in (3.21), but 
their sums are restricted to combinations of expansions which commute with each other. 
The first term involves integrals performed over the complete domain D and will be kept 
for the final result. The second term is reinserted on the left-hand side when the recursion 
formula (5.13) is iterated from n = 1 as in (5.6) to n = minjA'c -|- 1, — 1}. The recursion 
ends when n = Nc + 1 because this is the maximal number of regions which can be combined 
into a subset with commuting expansions, i.e. subsets containing all regions in Rc plus one 
region out of R^c each. In the last step of the recursion, when n = A'c + 1, the second term 
on the right-hand side of (5.13) is absent. 

Finally the third term on the right-hand side of (5.13) sums over pairs of regions 
with non-commuting expansions plus an arbitrary number of commuting expansions. The 
integrals in this term are performed over the convergence domain D^/^ of the non-commuting 
expansion applied last. Also these terms are kept unchanged for the final result, so no 
combinations of more than two non-commuting expansions appear in this formalism. 

The recursion results in the following expression for the original integral F: 

{Rc+l) {Rc+l) 

p = ^ pi^) _ ^ _p{^'i.a:2) + . . . — ( — 1)" ^ _p{a;'iv,0 
xeR {x[,x'2}cR {x[,...,x'„}cR 
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_|_ . . . _|_ ( — 1)^'= ^ ^{x',xi,...,xjvj 

x'£Rnc 

F2J / 

This is the generalized master identity for the expansion by regions with non-commuting 
expansions. The terms in the first two lines of (5.14) look familiar when compared to 
the ones from section 3 in (3.22). The first term is exactly the same: single expansions 
according to all regions, integrated over the complete integration domain, corresponding 
to the known recipe of the expansion by regions. 

Also the following terms in the first two lines of (5.14) are similar to the identity 
from section 3: Multiple expansions integrated again over the complete domain. Here only 
those combinations of expansions appear which commute with each other, but otherwise 
these terms are the same as in (3.22). Whether they are scaleless or relevant contributions 
depends on the same properties of the regularization and choice of regions as described in 
section 3. 

The terms in the last two lines of (5.14), however, are disturbing. Each of them is an 
integral performed over the convergence domain D^/^ of a non-commuting expansion. We 
would not expect such terms to appear from the expansion by regions. And we cannot be 
sure that the series expansions in these terms converge separately, because such restricted 
integrals are not covered by condition 4 from section 3. So we have to clarify why and 
under which circumstances these extra terms vanish. 

There is a general argument why these extra terms should cancel among themselves 
if we only have a small number of non-commuting regions such that the boundaries of 
their convergence domains can be varied independently. The complete integral F must 
be independent of the boundaries between the convergence domains D^. But all integrals 
in the large round brackets of (5.14) are only performed over the subdomain D^'^. If the 
boundary of D^'^ is varied infinitesimally without changing the boundaries of the other 
regions from Rnc, then the expression in the brackets and therefore the complete result 
changes — unless the sum of all the integrands within the brackets vanishes close to the 
boundary. And because the position of the boundary is quite arbitrary within some range, 
we may expect that the sum of the integrands vanishes all over D^'^ in order to make the 
complete result independent of the boundary. 

In the following we use a different, more formal way to show that the extra terms 
cancel. We impose an additional condition on the expansions: 

5. For every combination of two non-commuting expansions, there is a region from 
whose expansion does not further change the doubly expanded integrand: 

yx[,x'2 G Rnc, x[ / 4, 3xeRc: r(^2^4-^') = T^'A^<^ . (5.15) 
In every example with non-commuting regions which I have studied so far, this is indeed 
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the case. A few such examples are presented in sections 6 and 7 and in appendix B.2 of 
this paper. 

If condition 5 holds, then, for every pair x'^,X2 G Rnc, the extra terms in the round 
brackets of (5.14) can be grouped in pairs of two, with and without this particular region 
X € i?c involved: 

-(-1) 2^ (^[41 -^[4] ]+••• (^-^^^ 



{x'l,...,x',[}cRc\{x} 



With the equality of the expansions imposed in (5.15), these combinations of integrals 
vanish at the integrand level, as 

p(x2-<— p(x2-<— Xj) j-,{x2<—X-^,X,Xj^,...,X^) J-,{X2<—X-^^,X-^,...,X^) /_ ^ „x 

^[4] -^[4] ' ^[4] ' 

and therefore all extra terms in the round brackets of (5.14) cancel with each other and 
vanish. 

We conclude: If condition 5 holds, the master identity for the expansion by regions 
with non-commuting expansions reduces to 

(-Rc+l) (Rc+l) 

p = ^ pi'-^) _ ^ _p{^'i.4) + . . . - (-1)" ^ p{^'i:-y,j 

xeR {4,x^}cR {x[,...,x'„}cR 

+ ... + (-1)^^ _p(x',xi,...,xjvj ^ 

x'GRnc 

which reproduces the identity (3.22) in section 3 with the difference that the sums over sets 
of expansions are now restricted to such combinations which commute with each other. 

Note that, as in section 3, it is possible to restrict the derivation of this identity to the 
leading-order term in each expansion, yielding the leading-order approximation Fq (3.27) 
for the integral F with the same restrictions on the sums over regions as in (5.18). 

6 Example: Sudakov form factor 

A toy example for non-commuting expansions can be found in appendix B.2, where it 
is discussed in the context of divergence problems with the series expansions. Let us 
study here as a "real" loop integral with non-commuting expansions the one-loop vertex 
correction to the Sudakov form factor: a three-point function with one massive exchanged 
particle where one leg is off shell by a large amount, whereas the other two are on shell 
(and massless, for simplicity): 



F = jBkl , with / = /1/2/3 and 

T 1 1 

' {{k + Pi)T' (A:2 + 2pi • A;)"i ' 
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Pi + k 



m 



P2 + k 



Figure 3. Vertex correction to tire Sudakov form factor. 



(6.1) 



with Pi = P2 = and {pi — p2)^ = —2pi ■ p2 = —Q^ < 0. This integral is illustrated in 
figure 3. We consider the integral in the Sudakov limit 



(5^ > 



(6.2) 



so we want to expand it in powers of rr? jQ'^ . We are particularly interested in the case 
ni = = n-s = 1, but we will use the propagator powers, especially n\ and n2, as analytic 
regulators. 

In the Sudakov limit the integral (6.1) is characterized by the two light-like directions p\ 
and p2 ■ It is convenient to parametrize any momentum k through the light-cone coordinates 



and the perpendicular components such that 

A: = ^ -I- /c"*" ^ -I- fc_L , with p\^2 • fe_L = . 



k\ 



(6.3) 



(6.4) 



In terms of these light-cone coordinates, scalar products are given by 



k-l 

the propagators read 
h 
h 



(6.5) 



{k+k- -kl + Qk+)^^ 



{k+k- - k\ + Qk- 



{k^k~ — kj_ — m?) 



2\n3 



(6.6) 
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and the integration measure (2.3) factorizes into 



oo 

2e e7E r 



°' = ^n^y^ ydfc-dfc-. (6.7) 

— oo 

When looking for the relevant regions of this integral, let us start with the known hard 
and collinear types: 

• the hard region {h), characterized by ^ k± ^ Q, with the expansion 

tW/,, = /,,, r(^),3 = f;(!^ K)^ , (6.8) 

^ j! {k+k- - A:^ )n3+j 

converging absolutely within D/^ = ^k £ D : k'j_^ m^j, 

• the l-coUinear region (Ic), characterized by k^ ~ m?/Q, k~ ^ Q and k± ~ m, 
with the expansion 

^ '^'^ - '^'^ ' ^ - . ^ , iTT?^ (Qfe-)-+- ' ^'-'^ 

converging absolutely within 

Die = {A; G : 1^+1 < Q A A^ < Q|A;-| A A^ < m^} , 

• and the 2-collinear region (2c), characterized by k~^ ^ Q^k~ m? /Q and k^ ~ m, 
with the expansion 

ri2c)j _ (^i)ji2 (-fc+fc-)^'^ (fcj)^'^ (2,) _ 

^^1-^.^1.^7^;^^ (QA+)ni+.i. ' ^ ^^2,3-/2,3, (6.10) 

converging absolutely within 

D2c = {keD\Di^: |A;-| < Q A Ai < Q|A;+| A A?! < n?) , 

where D = is the complete integration domain. The domains D^, Die, have been 
chosen to cover as much of each expansion's convergence domain as possible, but in such 
a way that they are non-intersecting. In particular, in the domain with k'j_ < and 
k\/Q <^ \k^\ ^ Q, both collinear expansions T^^'^^ and T^^"^) converge. Here this domain 
is attributed to Die and therefore explicitly excluded from D2c-, but this is an arbitrary 
choice, and the boundary between Die and D2c can be chosen differently. 

For the convergence of rc^) we assume that |A;+A — kW ^ rm? can be obtained 
whenever fe^ S> rn?, through an appropriate bending of the integration contours away 
from the real axes where necessary to avoid cancellations between the terms k~^k~ and 
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Figure 4. Convergence domains of the regions for the Sudakov form factor in the case k\ < m} . 



k'j_. Such convergence issues have been checked expUcitly for the threshold expansion in 
appendix A. 2.1.'^ 

Before we search for further relevant regions, let us state why we do not need a soft 
region here. The soft scaling is /c"*" ~ A;~ ~ A;_l ~ m, and the corresponding expansion 
is equivalent to the double collinear expansion, T^^^''^^\ expanding both propagators Ii 
and I2 according to (6.9) and (6.10).^ So a soft, i.e. double collinear expansion is involved 
anyway, even if we do not add a soft region. Adding a soft region would not change the 
result within our formalism. As T^*) = T^^^'^'^) = T^i'^-^) = T^^^.^) = t'''^c,2c,s) ^ additional 
terms would cancel each other: 

(^ — l)"" ^pi^^^l'---''^") — ^(lc,s,xi,...,x„) _ j^{2c,s,xi,...,x„) _|_ jp(lc,2c,s,xi,...,x„)^ _ q |-g -^-^-^ 

for any set of other regions with Xi {Ic, 2c, s}. 

To determine which further regions we need, we take a look at which parts of the 
integration domain D = M.'^ are already covered by Dfi, Die and D2c- The domain with 
A;^ 3> is identical to D^. For /c^ < m^, the partitioning of the |/c+|-|/c~| -plane is 
schematically shown in figure 4. Two corners of this plane are not covered by Die and D2c- 
The domain where both and \k~\ are small, of the order k'j_/Q. And the opposite 
corner where both are large, of the order Q. We associate these two missing domains 
to convergence domains of additional regions and establish the corresponding expansions. 
To the list of regions we add 

^Based on the same reasoning, we could have included into Dh the domain with S> and 

fcx ^ "1^. But this only works as long as the term k^k" is present in the denominator of the third 
propagator. We also need a region where fc""" and k~ are both small (see below), expanding I3 by eliminating 
k'^k" from its denominator. A second expansion according to the hard region is then impossible within this 
domain where k\ < m^. So we choose to exclude this domain from Dh and rather establish an additional 
region especially dedicated to large fc^. 

*See the expanded integral (A. 47) in appendix A. 3.1 for details. 
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the Glauber region (g), characterized by A;+ ~ A; ~ m? /Q and ~ m, with the 
expansion 



r(.)/3 = f (-^^^^")^' , (6.12) 

converging absolutely within Dg = {k £ D : Q\k^\ < kj_ < rn^} , 

and the collinear-plane region (cp), characterized by /c"*" A;~ Q &iid /c^ ~ m, 
with the expansion 



rr{cp)r _ Y> (^l,2)j (fcj)^' 



i=o 

converging absolutely within Dcp = |A;GZ): |fc^|^QAA;^< m^}. 

While the Glauber region has some physical meaning and can yield relevant contribu- 
tions for other integrals or with other regularization schemes, the collinear-plane region 
is completely artificial and just needed to cover the integration domain with convergence 
domains. 

We now have five regions with non-intersecting convergence domains which fulfill 

D,,UDi,UD2cUDgUDcp = D. (6.14) 

Most of the corresponding expansions commute with each other, as one can easily verify, 
with one exception: The Glauber and collinear-plane expansions do not commute, their 
double expansions read 

Ticp^9)j= Y (^2)i34 (nsh. (-fc+fc-)J-i3^ {kly^^ 

j^^^^^O il' • ■■35^- {Qk+Y^+n2 (Qk~)n2+j34 (_A;2 _ m2)^3+i5 ' 

(„) jni),,, (n2),34 (n3),aa (-fc+fc-)^» (fcj)^'^^- (m^)^^ 

. ii!---j6! (QA;+)'^i+Ji2 (Q/fc-)"2+i34 (A;+A;-)"3+i56 • ' 

So we have a situation where the formalism for non-commuting expansions developed in 
section 5 is needed. The sets of commuting and non-commuting expansions are Rc = 
{h,lc,2c} and R^c = {9,cp}, respectively. Condition 5 (p. 28) is easily satisfied, there 
are always several regions whose expansions do not further change the doubly expanded 
integrands: 

rp{cp^g,lc) _ rp{cp^g,2c) _ rp{cp^g) 
rp{g^cp,h) _ j,[g4^cp,lc) _ j^{g^cp,2c) _ j^(g^cp) _ (6.16) 



-33- 



Therefore the extra terms in the identity (5.14) vanish: 

_ picp-^g) I p{cp-i-g,h) p{cp^g,lc) p{cp^g,2c) 
[cp] [cp] [cp] [cp] 

_ p{':p<-g,h,lc) _ j^{cp^g,h,2c) _ ^(cp-i-g,lc,2c) p,{cp<-g,h,lc,2c) 

[cp] [cp] [cp] [cp] 

= (_1 + 2 - 1) F^^-^) + (1 - 2 + 1) F^^^^''^) = (6.17) 

and 

_ rp(g^cp) rp{g^cp,h) p{g4-cp,lc) p(g^cp,2c) 

^[g] ^^[g] ^^[g] +^[9] 

_ rp{g^cp,h,lc) _ j^{g<-cp,h,2c) _ rp{g<-cp,lc,2c) rp{g^cp,h,lc,2c) 

^[g] ^[g] ^[9] ^^[g] 

= (-1 + 3-3 + 1)Fj5J^'^^ =0. (6.18) 

We are left with the terms of the identity (5.18), 

F = F^'^^ + F^^^^ + F^'^'^^ + F^^^ + F^^P^ - ^F^'^'^^^ + F^'^''^^^ + F^'^'^^ + F^'^'^^^ 

_j_ _p{lc,2c) _|_ p(lc,g) _|_ _p(lc,cp) _|_ p(2c,g) _|_ ^(2c,cp)^ 
_j_ p{h,lc,2c) _|_ p{h,lc,g) _|_ p{h,lc,cp) _|_ p{h,2c,g) _|_ p{h,2c,cp) _|_ p{lc,2c,g) _|_ ^(lc,2c,cp) 
_ ^^(^,lc,2c,g) _|_ ^(/x,lc,2c,cp)^ ^ (6.19) 

summing over aU combinations of regions where the expansions commute with each other. 
The evaluation of these contributions is described in appendix A. 3. The first subsection, 
appendix A.3.1, shows that the contributions F^^^ and F^'^^^ from the Glauber and collinear- 
plane regions are scaleless, and the same holds for all overlap contributions with multiple 
expansions. Omitting these scaleless contributions, the integral F can be expressed as 
follows: 

F = F^'^) + f(^^) + f(2<^) . (6.20) 

Remember that we have used analytic regularization, and we actually need it for having 
all scaleless contributions well-defined. 

If we had introduced a soft region and added F^'^^ to the result (6.20), according to the 
recipe in section 1 for the expansion by regions, then nothing would have changed because 
F(*) = F^^'^'^'^) = 0. 

We also note that even if the collinear-plane region did not yield scaleless integrals, 
all contributions involving the collinear-plane expansion in (6.19) would cancel each other. 
This is because adding the hard expansion does not change the collinear-plane expansion, 
pih,cp) ^ p{cp) (unless applied after r(^^)), which can be seen in the expressions (A. 51) in 
appendix A.3.1. So the contributions involving T^'^p^ can be grouped into pairs 

f^picp,xu...,x„) _ p(h,cp,xi,...,x„)^ ^ Q ^ ^.^j^ {xi,...,Xn}C {IC, 2c} , (6.21) 

which cancel each other. 
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The contributions remaining in the identity (6.20) read 
{n3)j I 2W / DA; 



pile) ^ f> M, [ m{-k^y 



j=0 



based on the expansions T^^^ (6.8), r(i=) (6.9) and r(2c) (6.10), but rewriting the expres- 
sions in a Lorentz-invariant way, for the colhnear regions using (A.41). The evaluation of the 
hard contribution F^^'^ is straightforward with Feynman parameters (see appendix A. 3. 2). 
The colhnear contributions are calculated in appendix A. 3. 3. The integrals yield 

^ ^2.^n^^-^.n,,, ~ /_m2y r(n3+j)r(ni23-2 + 6 + j) 

r(ni) r(n2) r(n3) ^ V i! r(4 - ^23 - 2e - i) 

X r(2 - 7113 -e-j) r(2 - n23 - e - j) , 

r(n0r(n2)r(n3)^"'^ ,t^V 

r(n2 + j) r(T2i - 712 - j) r(T2i3 - 2 + e - j) r(2 - Til - e + j) 

j!r(2-722-e-i) 

7?(2c) _ ^ E / 2N2-n23-^ ('02N-ni / _ !!L ] 

-r(r70r(772)r(7.3)^"'^ ^ V 

r(ni + j) r(722 - ni - j) r(7i23 - 2 + e - j) r(2 - 712 - e + j) 
j!r(2-7zi-e-i) 



(6.23) 



Two remarks on these results are in order. First, all terms in the series expansions in (6.23) 
are homogeneous functions of rv? and , where the leading-order terms scale as 

Fo(")oc(Q2)f-"i23^ 

Fo^'^) oc (7n2)i-"- ((fr^\ 

F(S^'^ OC (7^2) f -"23 (Q2)-ni _ 24) 

These scalings could have been obtained directly from the integrals (6.22), taking into ac- 
count in each region the scaling of the integration measure and of each propagator with 77i2 
and . Using such scaling arguments, the contributions from the Glauber and collinear- 
plane regions would scale as 

FjS'^) OC (7n2)f+i-"i23 (q2)-1^ 
F'^'"^ OC (77x2)1-1 (g2)l-ni,3 ^ (0.25) 
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if they were not scaleless. Each region's scahng has a unique dependence on the regular- 
ization parameters d and ni, 712,^3 by which it differs from the other regions. The overlap 
contributions originate from multiple expansions with different scalings, so we could have 
expected all overlap contributions to be scaleless even without calculating them. 

Second, the two collinear contributions F^^"^^ and F^'^'^^ in (6.23) are individually diver- 
gent when ni and n2 tend to integer numbers, in particular for the case ni = 712 = ns = 1 
we are interested in. The Gamma functions T{ni — n2 — j) and T{n2 — ni — j) , respectively, 
provide single poles in the variable (ni — n2) once the expansion order j is large enough 
(starting at j = for ni = 71,2). So here analytic regularization is not only needed to make 
scaleless contributions well-defined, but also to prevent contributing regions from yielding 
ill-defined integrals. 

Let us have a look at the particular case ni = n2 = = 1- We cannot simply set all 
propagator powers to integer values, but we have to understand this as a limit. We use 
the symmetric values ni = 1 + 6 and n2 = 1 — S and expand both collinear contributions 
in a Laurent expansion about (5 = up to the finite order 6^. This choice of the limit has 
the advantage that the overall dimension of the integral is not changed, but any other way 
to approach the point ni = 712 = 1 would yield the same result for the sum F^^'^^ + F^'^^\ 
though not for the individual terms. 

The hard contribution is not singular in the limit 6 — t- 0, so we can directly set 
77i = 772 = 773 = 1 here. According to (A. 54) and (A. 56) in appendix A. 3. 2, 



2 



1 



e^7Er(i + g)r^(-g) 
r(i - 2e) 

1 2 / 777^ 



m 



2F1 [ 2e, 1; 1 + e; g2 



2Li9 



771 



vr 
12 



.26) 



The original integral F (6.1) is finite in four dimensions (i.e. for e = 0) when 771 = 772 = 
773 = 1, because there are no ultraviolet divergences and the infrared (soft and collinear) 
divergences are absent due to the finite mass 777 of the exchanged particle. The hard 
contribution F'^^^ ^ however, involves infrared singularities (single and double poles in e) 
because the expansion, assuming large loop momenta, yields massless integrals. These 
singularities in (6.26) have to be cancelled by the contributions from the collinear regions. 

The collinear contributions are evaluated for n\ = \ ^ b and n2 = \ — b. The results 
obtained in (A. 62) and (A. 65) in appendix A. 3. 3 read 
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The 2-collinear contribution is obtained by replacing 5 — )■ —(5: 
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When the two colhnear contributions are added, the l/(5-singularities drop out. The total 
coUinear contribution is finite in the limit (5—7-0, i.e. in the case ni = n2 = n^ = 1: 
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The total collinear contribution exhibits the same single and double poles in e as the 
hard contribution F'^^'^ (6.26), just with opposite signs. So the sum of hard and collinear 
contributions is finite in the limit e — ?• 0: 



1 



1 



In^ 



+ In 



Q 



m 



m 



+ 0(e) . (6.30) 



This result agrees with a direct evaluation of the original integral F (6.1) for ni = 77-2 = 
723 = 1 via Feynman parameters without expanding in rr? jQ'^ . 

For general propagator powers ^1^2, 3 the following Mellin-Barnes representation of the 
full integral F (6.1) is found: 
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(6.31) 
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Figure 5. One- loop corrections to tlie forward-scattering amplitude with small momentum ex- 
change r. 

An asymptotic expansion for ^ is obtained by closing the Mellin-Barnes integral 
to the right and picking up the residues of the poles from the three T{. . . — z) functions in 
the numerator. If the difference (ni — 77-2) is non-integer, then each of these three Gamma 
functions provides a series of single poles, and the three contributions F^^\ fi'^'^) and F^"^^^ 
in (6.23) are exactly reproduced. Otherwise, especially in the case ni = 712 = 1, the poles 
of two of the Gamma functions coincide, providing a series of double poles. For this reason 
the contributions F^^^^ and F^"^^^ cannot be individually finite in the case ni = 712 = 1, but 
their finite sum, F^^'^^ + F^'^'^\ corresponds to the series of residues from the double poles. 
If we go one step further and set e = in (6.31) in addition to ni = 77-2 = ns = 1, then a 
series of triple poles arises which yields the finite result (6.30). 

The residues of the Mellin-Barnes integrand (6.31) are in one-to-one correspondence 
with the contributions from the (/i), (Ic) and (2c) regions. This correspondence links 
singularities arising in contributions from individual regions to coinciding poles of the 
Mellin-Barnes integrand. 

7 Example: forward scattering with small momentum exchange 

Let us look at a last example which permits us to study the dependence of individual 
contributions on different regularization schemes, in particular with and without analytic 
regularization. It also is an example where overlap contributions turn out to be relevant 
under certain circumstances. We evaluate a set of one-loop corrections to the forward- 
scattering amplitude of two light-like particles with a large centre-of-mass energy, but 
small momentum exchange between them. 

The two contributing diagrams are displayed in figure 5. For reasons to be seen later 
we symmetrize the loop integral under commutation of the momenta of the exchanged 
particles, /c -H- (r — k), before specifying general propagator powers ni_2,3,4- The loop 
integral reads 

F=^-[ ^ i I + ' ^ 



^ +T^—^-^T^\- (7-1) 



{{p2 + kYY' {{p^ + r-kY) 
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We are interested in the case ni = n2 = = = 1, but we will partly use the propagator 
powers as analytic regulators. 

All internal and external lines are massless. The on-shell conditions for the external 
particles, pf = {pi — r)^ = P2 = {P2 + ?')^ = 0, imply 2pi ■ r = —2p2 ■ r = r^. The centre- 
of-mass energy shall be much larger than the momentum exchange, = {pi + ^2)^ = 
2pi • p2 ^ |r^|. We use the light-cone coordinates specified in (6.3) and (6.4) such that 

= 2pi^2 ■ f/Q = '^r'^/Q and = r^r~ — f'j_. The only independent kinematical 
parameters are and r'j_ with 

Q2 > rl . (7.2) 
The dependent parameters are determined by 

...,1^.,1,0((|)!). 

In terms of the expansion parameter fJ_/Q'^, the exchanged momentum r thus has the 
scaling of the Glauber region. The loop integral is written using light-cone coordinates: 

CXD 



2 {k+k- - klY^ {{k+ - r+){k- - r-) - {k± - f^yy' 
1 1 




{{k+ + Q)k~ - kl)^ {{k+ -r+ - Q)k- - k+r~ -kl + 2r± • A;^) 
where relation (7.3) has been used to cancel certain terms in the denominators. 



(7.4) 



7.1 Regions for the forward-scattering integral 

One finds that the same regions as in the previous example of section 6 are relevant here, 
with m replaced by |rj_| for the scaling prescriptions. The convergence domains (x = 
h, Ic, 2c, g, cp) are similar, but slightly more involved. They are illustrated in figure 6, where 
the partitioning of the |A;"'"|-|A;~|-plane is shown for the cases |A;_l| ^ \ f±\ (left diagram) 
and |A;_l| < |rj_| (right diagram). 

In the present example we will study all expansions only to leading order in the ex- 
pansion parameter r^/Q^. In detail the regions we need are 

• the hard region (h), characterized by k^ k~ ^ k± ^ Q, with the (leading-order) 
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Figure 6. Convergence domains of the regions for the forward scattering with smaU momentum 
exchange. 
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where the all-order expansion converges absolutely within 
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(7.6) 



the 1-collinear region (Ic), characterized by k'^ ~ ^±/Q, k ~ Q k± ~ f±, 
with the expansion 



/ 1 1 

^7 + 



(A:+(A:- - Q) - fc2 )"3 (^+(^- ^ g) _ ^+^- _ ^.2 + • fc^)"'/ 



(QA;-)"4 (-Qfe-) 
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converging absolutely within 

i^ic = jfc G D : ^k^\ » \r±\ A < ^ A \k-\ » 

V < \r±\ A \k+\ « Q A \k~\ » ^^HM A |A;+| < |A;-|^ | , (7.8) 

the 2-collinear region (2c), characterized by A;+ ~ Q, A;~ ~ r'^/Q and k± ~ r^, 
with the expansion 

^i2c)j 1 1 1 



1 1 

+ 



(-QA:+)"3 (QA:+)"3 
1 



V {{k+ + Q)k- - kiy {{k+ - Q)k- - k+r^ -kl + 2r-l • k^^' J 

(7.9) 

converging absolutely within 

D2c=f^keD : ^k^\ » |f^| A \k-\ < ^ A l^+l » 

V (^k^\ < \f±\ A « Q A |A;+| » ^^^^M a > | , (7.10) 

• the Glauber region (gr), characterized by /c"*" ~ A;~ ~ ^±/Q and /c^ ~ rii with the 
expansion 



° ' 2 (_^2)n, (-(fc^-r-l)2)"^ 
/ 1 



\"4 



(7.11) 



\ {Qk~ - kl)"^ {-Qk~ -k\ + 2fi_ ■ ki 
converging absolutely within 

Dg = \^k(^D : l^k^\ » |f^| A < ^ A < 

y{\ki.\->\ri_\ A \k-\<^ ^ l^^l^^) 
v(M<|f.|A|.-|<^)l (7.12) 
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and the collinear-plane region {cp), characterized by /c+ ~ /c ~ Q ^^(^ k± ~ ri; 
with the expansion 



converging absolutely within 

Dep = {A: G L> : \k^\ < |fx| A > q} , (7.14) 



where D = M'^ is the complete integration domain and = ^r'^/Q is the leading-order 
approximation of according to (7.3). 

Note that without analytic regularization, i.e. setting ni = 77,2 = ^3 = n4 = 1, the 
Glauber contribution Fq^^ = jDkTQ^^I from (7.11) is only convergent because the integral 
has been symmetrized. Integrals over individual terms are divergent for \k~^\ — )■ 00 or 
\k~\ — )• 00 like dk^/k^, but the leading term in each of the round brackets of (7.11) is 
cancelled, making the integrals convergent at infinity like d/c^/(fc^)^. 

The convergence domains (x = h, Ic, 2c, g, cp) are non-intersecting and cover the 
complete integration domain, 

DhUDicUD2cUDgUDap = D. (7.15) 

It might look strange that the domains of the collinear and Glauber regions extend into 
the zone with |A:^| ^ |rj_| (cf. left diagram of figure 6). But this is necessary if we do not 
want to invent new artificial regions, because the hard convergence domain cannot cover 
the complete |/i;^|-|A;~ |-plane with \ k±\ ^> |ri|. The definitions of the convergence domains 
given above ensure that for every region x and every point k G Dx there is at least one 
large term left by T^^^ in the denominator against which the smaller terms according to 
the scaling of the region x can be expanded. 

The multiple expansions are determined from the scaling prescriptions for the regions 
given above; they are listed in appendix A.4.1. All expansions commute with each other 
with the exception of T^s) and T^'^p); 
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The expansion T^^^ does not alter the integrand when it is apphed after T^'^p\ because 
the Glauber domain Dg also contains a part with \k±\ ^> |rj_| (so possibly > Q) and 
\k~^\ ~ k\/Q or |A:~| ~ kj^/Q. Thus we cannot rely on |A:^| <C Q throughout Dg and 
have to leave {kr^ it Q) unexpanded in T^d-^'^v) (7.17). When, however, T^^-* is applied to 
the original integrand (7.4) or to the integrand expanded according to any other region 
except (cp), the term k~^k~ in the denominators is always accompanied by k'j_. So we can 
use |A;~''A;~| <^ k'j_, which always holds within Dg, in order to perform the expansion T^^^ 
as stated in (7.11). All other multiple expansions converge within the convergence domain 
of the last expansion, whatever expansions have been applied before. 

With the notations from section 5, the sets of regions with commuting and non- 
commuting expansions are Rc = {h,lc,2c} and R^c = {g,cp}, respectively, as in the 
previous example of section 6. Condition 5 from p. 28 holds because the doubly expanded 
integrand (7.16) is unchanged from further applications of T^^'^^ or T^'^'^\ and the expres- 
sion (7.17) is invariant under the hard expansion T^^^: 

rp{g4^Cp,h) _ rpig<-Cp) _ (7.18) 

So the identity (5.18) holds and has the same form as the identity (6.19) of the previous 
example. Here we restrict ourselves to the leading-order contributions: 

77 _ r^W I p(lc) p(2c) p(cp) / p(/i,lc) , p(/i,2c) r,{h,g) rp{h,cp) 



p(lc,2c) p(lc,3) p{lc,cp) rpi2c,g) r?(2c,cp) 

p(/i,lc,2c) p{h,lc,g) p(/i,lc,cp) rp{h,2c,g) rp(h,2c,cp) p{lc,2c,g) p{lc,2c,cp) 
p{h,lc,2c,g) _|_ p{h,lc,2c,cp)"^ 

Then, as explained at the end of section 3, we just have to ensure that these leading-order 
integrals are properly regularized (condition 3 of p. 14). And condition 4, concerning the 
convergence of the expansions outside their corresponding domains, is irrelevant. 

For the same reasons as in (6.21) we see that all contributions involving the collinear- 
plane expansion t'^^'p'^ cancel each other. 

Using scaling arguments, i.e. considering the scaling of the integration measure and the 
propagators in each region, we predict how the leading-order contributions of the individual 
regions depend on and (this will be confirmed by explicit evaluations later): 
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When the propagator powers Ui are used as analytic regulators, each region exhibits a 
unique dependence on and Q^, so we expect that all overlap contributions must be 
scaleless. Without analytic regulators, however, i.e. for ni = n2 = = = 1, the 
scalings read 
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(7.21) 

We notice that the hard and collinear-plane contributions, Fq'^^ and Fq^^\ are suppressed 
with respect to the "true" leading-order contributions Fq , ^o^'^^ ^^'^ ^o^^ ■ And the 
leading collinear and Glauber contributions share exactly the same dependence on 
and Q^. So there is no reason why the overlap contributions Fq^^''^^\ Fq^'^'^\ F^'^'^^ and 
p^c,2c,g) gj^Q^^^ scaleless; we must be prepared for relevant contributions from these 
multiple expansions. 

Altogether we expect that the leading-order result for the complete integral F (7.4), 
with or without analytic regulators (in any case for rij w 1 Vi and e ~ 0) is given by 

Fo = Fo^'^^ + F^^""^ + F^'^ - (Fq^'^'^^) + Fo^'^'^) + F^^""''^) + Fq^'^'^'^'^) , (7.22) 

where omitted contributions are either suppressed or scaleless, which is checked explicitly 
in appendix A. 4. 

7.2 Evaluation with analytic regulators 

The integrals contributing to (7.19) are evaluated in appendix A. 4.1 using the propagator 
powers ^1^2,3,4 as analytic regulators. It is shown there that the Glauber contribution Fq^\ 
the collinear-plane contribution F^^'^ and all overlap contributions are scaleless. For the 
non- vanishing contributions F^\ F^'^'^ and fI^^\ the scaling with r'j_ and predicted 
in (7.20) is confirmed. In the limit Ui — )■ 1 Vi, i.e. when all propagator powers tend to 1, and 
close to d = 4 dimensions (e 0), the hard contribution F^^'^ is suppressed with respect 
to the collinear contributions by one power of f'j_/Q'^. So the general expression (7.22) for 
the leading-order result is confirmed. When analytic regularization is employed, all but 
the first two terms are scaleless, and the leading-order result reads 

Fo = + F^^"^ . (7.23) 

These two collinear contributions are given in (A. 77) and (A. 79) of appendix A.4.1. They 
both involve a factor which is singular for ^3 = 77,4, but the limits ni — )• 1 and n2 — )• 1 are 
well-defined. Let us evaluate the contributions for ni = n2 = 1, 71.3 = 1 + ^3 and 714 = I-I-54: 

enBY{i + e + 53)T\-€-63) 



r(l + (^3)r(-2e-53-54) 
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e^7Er(l + e + <54)r2(-e-54) 



r(l + (54)r(-2e-53-(54) 



(7.24) 



We are interested in the case ni = n2 = = = 1, i.e. in the point 63 = 6^ = 0. The two 
colhnear contributions individually depend on the direction in the (53-54-plane by which 
we approach this point. One possible choice is the antisymmetric case of approaching the 
desired limit on the line 63 = —64^ = 6: 



e^7E r(i + g + ,5) r^{-e-6) 

r(i + 5)r(-2e) 

1 in //i^ Y e''^Er(l + e)r2(-e) 



2 f2Q2 \^f2 J r(-2e) 



r(i-5)r(-2e) 



2flQ^\rlJ r(-2e) 



. (7.25) 



In the limit 5—7-0 both collinear contributions are equal. Note that each collinear con- 
tribution is individually finite, which is due to a cancellation of 1/5 singularities between 
the two diagrams shown in figure 5, such that the factor (1 — e^^'^^) cancels the singularity 
from r(ib25). The total result with rii = n2 = = = 1 reads 

■ lime)L ..+ lime)l _^r/^^^e-^(l+^)r-(-e) 



(7.26) 

Another choice for the analytic regularization consists in taking the limit ^3 — )• first: 
The total result is then reproduced entirely by the 1-collinear contribution: 



Other choices are possible, but they all lead to the same total result. In fact, the analytic 
regularization parameters can be expressed through ^3^4 = A zt 5, isolating the singularity 
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in the variable 6. The total result remains invariant upon subtracting the 1/6 terms from 
each collinear contribution because these subtraction terms cancel each other: 
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(7.29) 



Here the expression in the square brackets has a well-defined finite limit for 5—7-0 and 
A —7- 0, independent of the way in which this limit is approached, which reproduces the 
known result: 
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(7.30) 



It is evident that lim„.^i Fq must be independent of the order in which the propagator 
powers ^^l,2,3,4 are sent to 1, because expression (7.23) is derived from the identity (7.19) 
which is valid for any regularization scheme (as long as the conditions of the formalism 
presented in sections 3 and 5 hold). 

The result obtained with the expansion by regions can be checked by evaluating the 
original integral F (7.1) with the help of a Mellin-Barnes representation: 
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(7.31) 



The asymptotic expansion for |r^| ^ is obtained by closing the Mellin-Barnes integral to 
the right and extracting the residues of the three Gamma functions in the last line of (7.31). 
The (possibly) leading terms originate from the residues at z = 0, z = 2 — 72123 ~ £ and 



2 — rei24 — £• Using 



^2 



'j_ and Q + r « Q , they reproduce the contributions 



from the hard (A. 68), l-coUinear (A. 77) and 2-collinear (A. 79) regions, respectively. 
7.3 Evaluation without analytic regulators 

Now we would like to see how the contributions from the previous section 7.2 change when 
we evaluate them without using analytic regulators. This means that we set all propagator 
powers to the fixed values ni = n2 = = = 1 before performing the integrals. 

The evaluation without analytic regularization of the integrals contributing to the 
complete expression (7.19) is described in appendix A. 4. 2. The scaling of the contributions 
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(A.91), F^^""^ (A.97), F^'^''^ (A.99) and F^^^ (A.lOl) with the parameters f]_ and 
as predicted in (7.21) is confirmed by exphcit calculation. This implies that the hard 
contribution Fq^^ is suppressed by one power of r'^/Q'^ with respect to the collinear and 
Glauber contributions and does not contribute to the leading-order result Fq. 

The collinear-plane contribution i^g'^^'' (A. 107) and all overlap contributions involving 
the collinear-plane region are found to be scaleless via dimensional regularization, even 
without analytic regulators.^ Also the overlap contributions involving the hard region 
yield scaleless integrals (A. 111). 

These results confirm the identity (7.22) for the complete leading-order result Fq, made 
up of the contributions from the 1-collinear, 2-collinear and Glauber regions including all 
their overlap contributions. The corresponding integrals are evaluated in (A.97), (A.99), 
(A.lOl), (A.103) and (A.106). They all yield identical results, 



77(lc) _ p(2c) _ _ p{lc,2c) _ p,{lc,g) _ p(2c,g) _ j^{lc,2c,g) 

,2 



rlQ^KflJ r(-2e) 



(7.32) 



This is very different from the contributions with analytic regularization determined in 
section 7.2. In particular, the terms (7.32) involve finite and relevant overlap contributions. 
The complete result, to leading order in r]_/Q'^, reads 

-TO — -Tq -1- -Tq -I- -Tq I-Tq -I- -t- ) "T -t^O 



(7.33) 



in ffi^Y e^TEr(l + e)r2(-e) 



f2Q2 yf2j r(-2e) 



The total result obtained without analytic regularization agrees with the result from the 
previous section, expressions (7.26), (7.28) or (7.30), where analytic regulators have been 
used. But the partitioning of the contributions among the regions is completely different. 
If we had omitted the overlap contributions from (7.33), the result would have been three 
times too large. 

At first sight it might seem to be just by chance that the cancellations between the 
various contributions in (7.33) leave exactly once the individual contribution. But this 
is a general feature whenever a certain set of regions (with commuting expansions) yields 
identical contributions, including all overlap contributions from these regions. Assume that 

p{x-^,...,x^) _ j Y)]^^^^^^ . . . J^^^")/ 

= Fo Vx;,...,x^ G = {xi,...,^^}, l<n<N<N. (7.34) 



subtle point concerning the regularization of poles pinching the integration contours at fc* = is 
treated in appendix A. 4. 2. 
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Then the total (leading-order) contribution from these regions, including the multiple ex- 
pansions, reads 

TV 



{x'^,...,x'^) 



n=l 



{x'^,...,x'„}cR 



Bk 



N 



1 



1=1 



I 



N 



1 



rid 

i=l 



Fn — Fn 



(7.35) 



where we have used /DA; (1—1) / = (1— l)i^ = = (1— l)Fo. The individual contribution Fq 
is left exactly once from all the cancellations. 

The Mellin-Barnes representation (7.31) is also valid without analytic regularization. 
For ?ii = 722 = n3 = = 1 it reads 



^2e ^e-y-E 



-2-t-z 



-2-e-z 



X r2(l + 2) r(2 + e + z) T{-z) r^{-l-e-z). 



(7.36) 



Closing the Mellin-Barnes integral to the right, the leading-order asymptotic expansion in 
IQ'^ is extracted from the residue at 2; = — 1 — e. The z-integral has a double pole there. 



and its residue yields the leading-order result (7.33), using 



and Q2 + ^2 ^ q2 



8 Conclusions 

Let us summarize the main statements of the general formalism introduced in section 3 
and generalized in section 5. 

The integral F = jBkl with integration domain D shall be expanded in some limit. 
We have defined a set R oi N regions, R = {xi, . . . ,xn}. Each region x is characterized 
by an expansion T^^^ = Ylj Fj^^ which converges absolutely when the integration variable 
is restricted to the corresponding domain D^- The conditions for the applicability of the 
formalism read: 



1 . The convergence domains are non-intersecting, H D^' 
complete integration domain, (J^ei? -^i' ~ ^• 



Vx 7^ x', and cover the 



2a. All expansions corresponding to regions within some subset Rc C R commute with 
each other and with expansions of any other region in R: 

rp{x)rp(x') ^ ri.(x')j,{x) £ , x' R . 

We write Rc = {xi, . . . ,XAr^} with Q < Nc < N and i?nc = R\Rc = {xNc+i^ ■ ■ ■ 
3. The original integral F and all integrals over expanded terms. 



j-,{x'^,x'2, 
jl,j2,- 



Jl 32 



are well-defined within the chosen regularization scheme, even if these integrals are 
not restricted to any convergence domain Dx- 
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4. The series expansions 

M,x'^,...) _ V- 

converge absolutely, although the expanded terms have been integrated over the 
complete integration domain D. 

5. For every combination of two non-commuting expansions, there is a region from Rc 
whose expansion does not further change the doubly expanded integrand: 

Vx;,4 G Rnc, x[ ^x'2,3xgRc: t^^)t^^'2)t^<) = r(4)rK) . 

Note that this condition requires Nc'>'i. 

If all these conditions hold, then the original integral is reproduced exactly, i.e. to all orders 
in the expansion, through the master identity (5.18): 

{Rc+l) {Rc+l) 
_p = ^ _ ^ _|_ _ _ _ _ ^ ^K>-.0 

x&R {a;;,x^}Ci? {x'-^,...,x'^}(ZR 

+ ... + (-1)^^ ^(x',xi,...,xjvc) ^ (8_1) 

where the summations with superscript {R^ + 'i) are defined in (5.8). They are restricted to 
such subsets of regions with at most one region from R^c- This means that the identity (8.1) 
sums over exactly those combinations of regions whose expansions commute with each 
other. 

The first term on the right-hand side of (8.1) corresponds to the known recipe of 
the expansion by regions presented in section 1. The following terms represent overlap 
contributions originating from multiple expansions. These overlap contributions may be 
relevant under certain circumstances, see the example from section 7 evaluated without 
analytic regularization. In many cases, however, the usual recipe of the expansion by 
regions is recovered, especially when dimensional and analytic regularization are used. If 
the contributions F^^-* with single expansions yield homogeneous functions of the expansion 
parameter in each order of the expansion and if this dependence on the expansion parameter 
is unique for each region, then we expect all overlap contributions to be scaleless, and only 
the usual contributions from each region survive. 

This means that in situations where the expansion by regions has normally been 
employed the known recipe [1-4] remains correct. The original authors of the method 
have always stressed how important the homogeneity of the contributions is. Where non- 
homogeneous contributions from individual regions appear, the formalism presented in this 
paper may still be applied if its conditions hold, and overlap contributions might arise. But 
it is usually preferable to change the choice of regions or regulators in order to obtain only 
homogeneous contributions and get rid of overlap contributions. In any case, one can al- 
ways check whether overlap contributions are relevant by simply evaluating them to see if 
they vanish. 
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The original recipe of the expansion by regions [1-4] can be understood such that 
any scaleless integral must be set to zero, whether it is well-defined through regularization 
or not. However, this is not the approach of the formalism presented here. Throughout 
the whole paper only such scaleless integrals are dropped which are mathematically well- 
defined through some regularization (and analytic continuation) and which can explicitly 
be shown to vanish (see appendix A. 1.2 and other appendices). Wherever an integral is 
ill-defined, a suitable regularization is added. Vanishing scaleless integrals thus are not 
a requirement for this formalism, but simply a well-understood property of dimensional 
regularization and analytic continuation which is used here. 

The leading-order asymptotic expansion of the integral F is obtained by replacing 
each series expansion in (8.1) with its leading-order term. Some of the leading-order terms 
obtained in this way may be suppressed with respect to others and must be dropped from 
the leading-order result. For the leading-order approximation condition 4 above does not 
apply, and condition 3 only needs to hold for the leading-order integrals. 

The examples presented in this paper illustrate the application of the formalism to var- 
ious integrals. This implies choosing the regions with corresponding expansions, covering 
the integration domain with the convergence domains and checking the conditions listed 
above. Most example integrals have alternatively been evaluated and expanded with the 
help of a Mellin-Barnes representation. For these simple one-loop examples, the asymptotic 
expansion via Mellin-Barnes integrals can be easier than the application of the expansion 
by regions. My experience with more elaborate multi-loop integrals depending on several 
small and large parameters, however, shows that the asymptotic expansion via multiple 
Mellin-Barnes integrals is often more difficult than the expansion at the integrand level 
within the strategy of regions. 

What is the practical use of the formalism presented in this paper, apart from the gen- 
eralization of the expansion by regions to cases where overlap contributions are needed? In 
principle, it should be possible also for more complicated integrals to prove their asymp- 
totic expansion by checking the conditions listed above. Most of these conditions can 
indeed be checked in a straightforward way once the list of regions, expansions and con- 
vergence domains is established. On the other hand, the first obstacle in getting there is 
finding all relevant regions such that the convergence domains of their expansions cover 
the complete integration domain. And this step might be quite tedious for complicated 
multi-loop integrals. It can be easier to evaluate successfully the contributions from the 
regions than to prove, using this formalism, that the evaluation is correct and complete. 
In particular, the formalism may need extra regions with scaleless contributions in order 
to cover the integration domain with convergence domains or to ensure the commutation 
of the expansions. 

Even with the formalism at hand, we will certainly not always want to prove its 
applicability, if only for practical limitations. But I hope this paper will convince readers 
that the expansion by regions is a well-founded method of asymptotic expansion and that 
its applicability can at least in principle be proven case by case. The standard recipe of the 
expansion by regions remains valid if overlap contributions are absent. So the formalism 
and its conditions provide hints on the proper choice of regions and regularization schemes 
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which ensure the homogeneous and unique dependence of the individual contributions on 
the expansion parameter. Finahy, the present paper shows that the expansion by regions, in 
its generahzed form, can be apphed to very different kinds of integrals, even such with finite 
boundaries or scaleful regularization parameters where the contributions do not exhibit a 
homogeneous dependence on the expansion parameter. Appendix B illustrates such an 
example and points out possible convergence problems arising there. 

Having obtained such knowledge on the foundation and generalization of the expansion 
by regions, this method is ready to be employed trustfully to various kinds of asymptotic 
expansions. 
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A Technical details of the evaluations 
A.l Large-momentum expansion 

This appendix presents details of the evaluation of the off-shell large-momentum expansion 
in section 2. 

A. 1.1 Large-momentum expansion: soft contributions 

The soft-region integral 



in (2.9) can be treated as an analytic function of the summation indices ji and j2- We 
first assume Reji^2 < and later analytically continue the result to positive values of ji^2- 
Using alpha parameters 




(A.l) 




oo 



(A.2) 







the integral (A.l) is written as 



g-j7r(n2-ii2)/2 



oo 



r(n2)r(-ji)r(-j2) 



doi da2 das 



"2-1 ^~J2-l joi(-m2+j0) 

1 2 3 
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with the shorthand notation (2.6). The loop momentum is integrated via 



d^ 

d/2 



ITT' 



P 



a 



+ iO 



yielding the alpha-parameter representation 



(A.4) 



-i7r(n2-ii2+2-e)/2 



r(n2)r(-ii)r(-i2) 



daiaf-^e^"i(-™'+*°) / daa «; 



a 



12 



X / das "'^ ^ 



exp 



(A.5) 



of the soft-region integral. After substituting 03 = t^/^ and evaluating the t-integral by 
reversing (A. 2), the a2-integral is solved using 



da a''-' (A + Ba) 







r(p) 



(A.6) 



followed by the ai-integral via reversing (A. 2). The result reads 

^2eg.7E ^-^n{n2-h-j2/2) (^^2^2-n2-e+n+j2/2 ^_^2 ^ .q)J2/2 

r ( - f ) r (2 - 6 + ji + f ) r (n2 - 2 + e - ii - f 



2r(n2)r(-j2)r 2-e + 



• (A.7) 

With the help of the doubling formula for the Gamma function (see e.g. [14]) we rewrite 

r(-f) _2i+^-^0F 



r(-j2) 



(A.^ 



Now we can safely perform the analytic continuation to ji^2 ^ 0. The factor ^{^-^p-^ in the 
denominator makes the result vanish for odd 72, as it should due to the Lorentz structure 
of tadpole tensor integrals. For even j2 we can write 



1 



(-1)^2/2 _ (-1)JV2 



J2\ 



r(i) ^2;,2/2 ^ (^]! 



where the relation 



r( ■\ ( r(a)r(i - g) • r(Q) 

r(a - j) = [.-ly -f^Ti — ^ = (-1)-^ 



r(l -a + i) ^ {I -a), 
for integer j has been used. Inserting the result into (2.9) we obtain, for even j2, 

,2\ ii+i2/2 



(A.9) 



(A.IO) 



m 



r(2-e + ji + f )r(n2-2 + e- ji-^ 



(-l)i2/2r(^^+^-^2) 



r(ni)r(n2) 



ii!(f )!r(2-6 + f 



(A.ll) 
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and = for odd 72. Rewriting the sum over ?'2 as a sum over j = ji + ^ and 

•?2 = according to (2.13), we calculate 



^ j^!(i-j0!r(2-6 + i^) 



j dxx"i+-'-i (1 - x)!-"!-^--'' (^^^^ (- 



i!r(2-ni-e-i) . 

J2=0 



-X 



1^2 



(1-x)^ 

r(ni+j)r(2-ni -e) 



j!r(2-ni-e- j)r(2-e + i) 
and obtain the result (2.14) for Fj'^^ . 

A. 1.2 Large-momentum expansion: overlap contributions 

This appendix is dedicated to the extraction of ultraviolet and infrared 1/e poles from 
the scaleless overlap contributions (2.37) in the particular case ni = 1, n2 = 2 where the 
complete integral F is finite. We have to evaluate 

Ji!j2! (p2)i+^i2 J {k^Y+i-n ' 



i,ji,h=0 iji,j2=0 



which only consists of vanishing scaleless integrals, but involves contributions of the form 
(1/euv ~ I/^ir) which we want to extract. The loop integral can be solved in analogy 
to the previous appendix A. 1.1 by assuming Rej2 < before analytically continuing the 
result to j2 > 0: 



00 

J (fc2 + io)- r(n)r(-j2)y 1212 J 





00 00 

2 P 



e^^E / doi a^-3+^ / da2 a^^'" exp 



r(n)r(-i2) 







ian \- iO 



= /i^^ e^^^ (n+2-,2/2-e)/2 (_^2 ^ .q),2/2 V ^ / „"-3-,2/2+. ^ .^3^ 

2r(n)r(-j2) J 


The prefactor is rewritten with the help of (A. 8) and (A. 9), vanishing for odd values of j2, 
as it should. The ai-integral is again scaleless. For extracting the singularities at ai — )• 
and ai — 7- oo we have to separate the integration at some intermediate scale A. Consider 
the following integral for integer j: 

oo A oo 

daQ^-i+^+ daa^-^+'= =0. (A.14) 

Ree>— j Ree<— J 
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The first term is ultraviolet-divergent for Re e < —j while the second term is infrared- 
divergent for Ree > —j. The two terms are combined using analytic continuation, and the 
integral vanishes. We are interested in the separate 1/e singularities, which only occur for 
j = 0: 

oo 

/ daa^-^^' = djo(— -] . (A.15) 

J ' V^uv eiR/ 



Plugging (A.15) into (A. 13) we obtain 

= (— - -) (py"' — ^ (A^i6) 

{k^r Veuv eiR/ r(n)l4^l! 



for even j2 and zero otherwise. Inserting this into (A. 12) with n = 2-\-i — ji, eliminating ji 
with the Kronecker delta in (A.16) and replacing j2 = yields 

^ -vA-^-^Ji-M) ^ ^^^ ^ }-Jw^m^^- ^"--''^ 

The sum over J2 can be evaluated in the following way: 

eoJ^!(^-j^)!(l+J^)! ^-oil^U;^ ^ r(2 + j^ + 5) 



1 



J2=0 



(l-x)' 



= r(i + <5) ^ 1 ^ . 

5^0 r(i - i + 5) r(2 + i + 5) r(i - *) (i + i)! ' 

For fixed z > 1 all terms in (A. 17) cancel each other in the sum over J2, resulting in 
l/r(l — i) = 0. The only remaining contribution is from i = 0, and we obtain 

F(h,s) = J_ f J LV (A.18) 



Combining the complete result F = F^^'^ pi^) — pih,s) (2.35), the poles in F^'*'*^ cancel the 
corresponding ones in (2.17): 1/euv cancels the ultraviolet pole in Fq*^ and l/ejR cancels 
the infrared pole in F^'^^ . 

A. 2 Threshold expansion 

Details of the threshold expansion in section 4 are dealt with in this appendix. 
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A. 2.1 Threshold expansion: check of convergence 

Let us check that the expansions (4.5)-(4.7) of the threshold expansion in section 4 really 
converge within their corresponding domains: 

• In the hard domain Dh = {k £ D : \k\ ^ \p \ V j/col ^ |p Iji we either have \k\ ^> \p\, 
then |p • k\ < \p\ \k\ <^ k^. If this is not the case, then there is \k\ < \p\ <C I/cqIj 
which implies \p- k\ ^ |pp ^ ko^ol- So \p- k\ is small compared to either P or | (70^0 1 
in (4.5), and these two terms are present in the denominators even if rC') is applied 
to T^P^Ii^2 where /cq has been eliminated from there before, cf. T^^''P^Ii^2 in (4-9). 

It remains to be shown that \p ■ k\ is also small compared to the combination of 
terms | A:q — fc^ it go^ol) which could be spoiled by a cancellation among the individual 
terms. Keeping track of the infinitesimal imaginary part in the propagators expanded 
withT^, we look for zeros of the denominators (fcg — /c^±gofco + ^0) and (/cq — fc^ + iO) 
in the complex fco-plane where the integration is performed along the real axis. These 
zeros exclusively lie in the upper left and lower right quadrants of the complex plane. 
For \k\ ^ \p\ zeros on different sides of the integration contour are well separated from 
each other. So we can bypass all zeros by bending the contour of the A^o-integration 
away from the real axis near the zeros. For |A;| < \p \ the zeros may pinch the contour 
near ko = 0. But then Dh requires \ko\ ^ where only two zeros, far apart from 
each other, are found. The same argument also works for the case when T^^^ is applied 
to r(p)/i 2. Therefore we can always choose the integration contour within Dh such 
that \p- k\ < - P±qoko\ in (4.5) and \p- k\ < |-P ibgo^ol for T^^^Phi^2 in (4.9). 

• In the soft domain Ds = {^k ^ D : \k\ < |A;o| ^ \p\]-, we have k"^ ^ k^ < \p\ |/co| ^ 
Igo^ol and \p- k\ < \p \ \k\ <^ If^o^ol) so each of the factors in the numerator of (4.6) is 
small compared to the denominator. 

• In the potential domain Dp = \^k ^ D : \ko\ '^\k\ <\p\^ , there is k^ <^ k'^. For the 
expansion T^p^Ii^2 in (4.7), we additionally need fcg ^ | — P it qoko — 2p- k\. Again 
the two zeros of (— fc^ it qoko — 2p- k + iO) exclusively lie in the upper left and lower 
right quadrants of the complex ko-plane. The /co-contour is pinched at ko = in the 
limit I A; I — >• 0, but as Dp requires |A;o| <C \k\, we can discard this case when checking 
the convergence of r(p). Otherwise the two zeros only pinch the fco-contour when 
A;^ + 2p • A; — 7- 0, which, for finite \k\, requires a certain angular correlation between 
k and p. So the pinching can be avoided if the contour of the angular fc-integration 
is bent into its complex plane in order to bypass the zeros of 1? + 2p ■ k. The 
unpinched /cg-contour may bypass the zeros in the fco-plane by bending away from 
the real axis. If T^^^ is applied to T^^^Ii^2, then the feo-contour is only pinched in the 
limit |fe| — ?• 0. So the integration contours within Dp can always be chosen such that 
\k^\ <^\-P± qoko - 2p- k\ in (4.7) and <^ \-P ± qoko\ for r(^'P)li,2 in (4.9). 

A. 2. 2 Threshold expansion: hard contributions 

In the hard-region integrals of F^^^ (4.10), the /co-integration converges for positive 
ni, 712, 71.3. We do not need the analytic regularization here and set ni = 712 = = 1 
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from the start: 



31,32=0 ■ 



{kl -k^+ qoko + iOy+n (A;2 - k^ - qoko + iO^+J^ {k^ - k^ + iO) 

(A.19) 



The first two propagators are combined with Feynman parameters, adding the third prop- 
agator afterwards: 

_f^^e^ ^ (2 + ji2)! ^rl^A.,^h(^_^\j2„l+jl2 



^TT^/^ jjj2! J 



ii,i2=o 





oo _ 

(2p- A;yi2 



X /dfeo /d-^-^fc ^,^ ,^ . -3^. (A.20) 



-oo 



The /c-integral is factorized into a radial |fe|-integration and an angular integration. For 
the latter we use 

j dQd-l = y d cos e (1 - COs2 0)('^-4)/2 J dO^_2 , J dJ],„2 = ^(O^ ' ^^'^^^ 

to distinguish the direction of the external momentum p, writing p-k= \p\\k\ cos 9. The 
angular integral vanishes for odd ji2 , and for even ji2 the fe- integral yields 



1 



(2|p 1)^12 -^-^ [ d C0S2 e (C0S2 0)Ol2-l)/2 (1 _ pog! 

r(i -e) y 



2 /i\-e 







r(i±iii')rf^ + e: 



(/c2 - A:2 - (2x - l)ygoA;o - ^O) 



3+jl2 



The /cQ-iiitegral over the last factor in (A. 22) is performed with a shift ko ^ k^ = ko + 
{x - i) yqo: 



oo 2 



dk'n 



2 



/2wl/2 



° 9 . 1.2 „ o\(3+ii2)/2-H 

1/2 /'/'^_ 1^2,2 2 



((-l-zO)fe^,V(x-i)%2g2^ 



l-i0)"i/2 ^ ^A.23) 
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By introducing j = J12/2, the summation over ji,j2 with even ji2 is rewritten as 



1 2J 



X 



b r 2 1 ii=o 

^ V ' 

= [x+(l-x)]2j=l 

The X- and y-integrations are easily solved. By substituting 

(2j)! = ^r(i+j)i! (A.25) 
the result F^'^) in (4.14) is obtained. 



A. 2. 3 Threshold expansion: potential contributions 

The potential-region integrals of -F^^^ (4.10) diverge at — t- ±00 if J123 > 0. For evaluating 
the potential contribution to all orders in ji , j2 , js we therefore need to keep general complex 
propagator powers ni and n2 as analytic regulators. The first two propagators of F^P^ (4.10) 
are combined with Feynman parameters: 



j,ip) ^ ^-^"""^ V r(ni2 + Ji2)(n3),3 ^n.^n-i n - x) 

^^^72 r(ni)r(n2) . ^ ^ ji!i2!j3! J ^ ' 

_ 00 

(p)m+Ts_J^ ° (fc2 + 2p • fc + (1 - 2x)qoko - zO)"^'+^'^' 
Partitioning the fco-integral into the two intervals (— oo,0) and (0, 00), we obtain 

dko ^ 

^ (^{k^ + 2p-k-iO) + {l-2x-iO)qokoj 

00 

+ /d(-A:o) ^ 

i ({k^ + 2p-k-iO) + {2x-l-iO)qo{-ko) 



"2+^2-1 



(2ji23)! r(ni2 - 1 - jl2 - 2j3) 1 1 _ 

123 

. (A.27) 



r(ni2 + ju) {k2 j^2p-k- i0)"i2-i-ii2~2i3 gi+^ii 

1 1 



(1 - 2a; - i0)i+2ii23 (1 _ 2x + i0)i+2ii23 
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Under the x-integral, this term extracts the residue at x = i: 



1 

J ^ ' \{l-2x-i 



-2m Res 



-i0)i+2ii23 (1 _ 2x + i0)i+2ii23 



4J123 (2ji23)! \dtj 
with t = X — \. The A;- integration 



(1 - 2x)i+2ii23 



t=0 



(A.28) 



/ 



{k^ + 2p-k- iO)"l2-l-Jl2-2j3 (/i;2)n3+. 



J3 



(A.29) 



is obviously scaleless if ni = n2 = 1 and any of the ji,j2,j3 is larger than zero (cf. [1]). But 
exactly in this case the Gamma function in the numerator of (A. 27) would be ill-defined, 
so we still need analytic regulators and continue the calculation with general ^1,722,723. 
Combining the two denominators with Feynman parameters, the /c-integral yields 



r(ni23 - 1 - jus) 



r(ni2-l-ji2-2j3)r(n3+i3) 



vr('^-i)/2r(ni23-i + 6-ii23) 



yrii2-2-ii2-2i3 ^-j^ _ y-jns+js-l 

d'^-^k 



(y{k + upY ~ ip'v'^ — iO) 

(-^•2 _ ,^0)|-"123-e+il23 

1 

X 



"123 — 1— il23 



r(ni2 - 1 - jl2 - 2j3) r(n3 + J3 

y dyy3-"i2-2n3-2e+ii2 _ 


7r('^-i)/2 r(ni23 - I + e - J123) r(4 - - 2m - 2e + jia) 



"3+i3-l 



r(ni2 - 1 - ii2 - 2^3) r(4 - ni23 - 2e + J123) 



X (p^ - iO) 2 



-".123— e+jl23 



(A.30) 



This cancels the singularity from (A. 27) when [jiyi — 1 — jvi — 2j3) is zero or a negative 
integer. The complete contribution 



^2£ g£7E 



-J7rni23 



TT 



'(p2 -i0)"i23-5/2+e r(ni)r(n2) 



jl,j2d3=0 

a \ 2ji23 



(713)^3 r(7Zi23 - 1 + e - J123) r(4 - ni2 - 2n3 - 2e + ji2) 



ji! J2! j3!r(4 - ni23 - 2e + ji23) 



ni+ji-l ^1 _ ^^n2+j2-l 



t=0 



(A.31) 
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is regularized dimensionally. We can switch off the analytic regularization at this point 
and set ni = n2 = = 1: 



il J2 J3=0 

2jl23 

(A.32) 



t=0 



Now it is obvious that the potential region gets no contributions from ji23 > 0: The 
maximal power of the polynomial in t in the last line of (A.32) is (ji +J2)) so this polynomial 
vanishes upon the 2(ji + j2 + J3) derivatives with respect to t. The only non- vanishing 
contribution originates from ji = j2 = Js = and is reported in (4.14). 

The leading-order potential contribution (from ji = j2 = ja = 0) is well-defined 
without analytic regularization. So for checking the result we can set ni = n2 = = 1 
from the start. We solve the /cQ-integration by closing its contour in the imaginary infinity 
and picking up one of the two poles: 



pip) _ 



_ 00 

2ege7E fd'^-^k f dkf) 



ITT' 



'^1'^ J P J (qoko-P -2p-k + iO)(qoko + P + 2p-k-i0) 

—00 

1 /i^^e^'^E r dd-ij^ 



qo 7r'^/2-i J {k^ + 2p- k-iO)k^ 

(A.33) 



1 f Y e'^^ V^r(He)r(-2e) 



y/q^p^-iO) W-iOj r(l-2e) 

using (A. 30) with ni = n2 = = 1 and ji = j2 = js = 0. This agrees with the result 
obtained from (A.32). 

A. 3 Sudakov form factor 

This appendix sketches the evaluation of the contributions to the vertex correction in the 
Sudakov limit from section 6. 

A. 3.1 Sudakov form factor: scaleless contributions 

First we have a look at the contributions from the Glauber and collinear-plane regions 
and at the overlap contributions which all turn out to be scaleless within dimensional and 
analytic regularization. 

Glauber contribution. According to the expansion 

T^9) (6.12) the Glauber contribution 

reads 

^2.g.7E ~ (ni),,(n2),,(n3),3 , ,,n23 /■ ^'^"'^^ 



pig) = V'Hjji Vii'2)j2 l'^3;i3 ( -|Y 

■71 I I'o I 1.3 I 



2iT^d/2 Ji!j2!j3! y (-A?2 - m2 + i0)"3+i3 

CO 00 

^ 7 dfc+(fc+)^-^^3 r dk~{k~y^^^ 

J (-B + Qk+ + zO)"i+Ji J (-P + Qk~ + i0)"2+i2 ■ 



-00 
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The A:^-integrals can be rewritten via = Qk^ — k\: 



00 

/ 



00 



dk^ (A;±)J123 1 f ^ ^ (t± + fc2 )il23 



—00 

These t^-integrals are scaleless, however. Within analytic regularization, we may assume 
that Re ni > 1 + j23 and Re 712 > 1 + jrs (separately for every term under the sum 
^ ji j2 j'i ) ' ^^i'^'^ required for the convergence of the t^-integrals and also allows to 
close the integration contour in the upper half of the complex t^-plane. The only pole at 
= — iO and the branch cut slightly below the negative real axis are then located outside 
of the closed contour. So both t^-integrals and therefore the fe^-integrals vanish. The 
A;j_-integral is regularized dimensionally. We get 

F(5)=0. (A.36) 

CoUinear-plane contribution. From t'^^p^ (6.13) we get 

jl,...,j4=0 

00 

X / d^""^^" (A 371 



Ficp) ^ y (m).. in2h (n3).34 (^2),. f ^.-2^ (^2 ),..3 



-00 



{k^{k- + Q) + iO)"'^^' ((fc+ + Q)k~ + iO) (A;+A:- + iO)"3+i34 
Here the /i;_|_-integral is scaleless: 



y A^~''k^ {k\r^^ = ^1^-^ I d\k^\ \k^\'+^^^'^-^^ = 0. (A.38) 



In the I /c_|_ I -integral the singularities at |A:_|_| — )• and \k±\ — )• 00 are regularized dimension- 
ally (by e) and cancel each other. The /c"'"-/c~-integral is regularized analytically. It yields 
the contribution 

2 e-^""^^3 (_i)ji234 r(ni23 + ji234 - 1) r(l - ni3 - ^134) r(l - n23 - J234) 

(Q2)ni23+Ji234-1 r(ni+ji)r(n2+j2)r(2-ni23-jl234) ' ^' ' 

which is non-zero and well-defined for ni^2,3 ^- Altogether we have 

= . (A.40) 

Hard— coUinear overlap contributions. For a Lorentz-invariant integration, the 1-col- 
linear expansion T(i^) (6.9) can be rewritten as 

r(ic)r ^ V K),.2 {-k^k-)n ikl)^ - (n2), i-ky 
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where terms in the numerator with the same j = ji + j2 have been combined via the 
binomial theorem. Putting (A. 41) together with T^^^ (6.8), the overlap contribution reads 



J^{h,lc) ^ (?^2)ji (r^3)i2 ^ ^^2y2 



ii j2=o 



DA: 



f A 421 

{k^ + 2pi ■ A;)"i {2p2 ■ /c)"2+Ji {k'^)^s+j2-ji ■ ^ ■ ' 

We evaluate this loop integral for general propagator powers ni^2,3 with the help of alpha 
parameters (A. 2): 

Dk 



I 



713 



3 ^-iwni/2 



doj a 



DA; exp 



i(ai3A;2 + 2A:- (aipi + a2P2)) ■ (A.43) 



Performing the loop integration via (A. 4), the Q2-integral through reversing (A. 2) and then 
the ai-integral using (A. 6), a scaleless aa-integral remains: 



2e „e7E 



H e 



g-i7r(ni3+2n2+2-.)/2 ^{n^ - ^2) r(2 - Tli - e) 



dQ3 Q^^^-=^+^ 



. (A.44) 



(Q2)n2 r(ni)r(n3)r(2-n2-e) 



Plugging this result into (A. 42) with n2 — )■ n2 + ji and 723 — )■ 723 + j2 — ji, we obtain 

f(^'^'=)=0. (A.45) 

The two collinear regions are symmetric to each other via exchanging pi ^ p2 and rii ^ n2- 
Therefore we also have 



(A.46) 



Collinear collinear overlap contribution. Rewriting both collinear expansions 
T^^*^) (6.9) and T^^^^^ (6.10) via (A. 41), the collinear-collinear overlap contribution reads 

DA; (A;2)Ji2 



^{lc,2c) ^ V- (^i)ji (^2)i2 f 

,4^, Ji!j2! ^ ' J {2pi-kY^+n 



jl J2=0 
^ (m)ii (n2)i2 

ii,i2=o 

X 



Ji!j2! 



J12 

1=0 

DA; 



J12 



[m 



{2p2 ■ A;)"2+i2 (/j2 _ ^2-)n3 
2v'i2-« 



fA 471 

{2pi ■ A;)"i+ii (2p2 • A;)"2+J2 (^^2 _ rn?Y^-^ ' v • ; 

where (A;2)Ji2 = ((A;2 

— m2)+m2)"'^^ in the numerator is expanded according to the binomial 
theorem. Again using alpha parameters (A. 2), 



/(2pi-A;)"i 



DA; 



3 — i7rni/2 



(2p2 • A;)"2 (A;2 - m2)"3 111 T{ni) 



n 



X y"DA; exp|^i(^a3A;^ + 2A; • (aipi + 02^2)) • (A.48) 
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Performing the integrations over /c, ai and as, a scaleless a2-iiitegral remains, 

OD 

-«7r(3ni+n2+2n3)/2 „ _ o i r 

^2. g.7E ; / da2 a"^-"i-^ = , (A.49) 

^ (Q2)ni (^2)ni3-2+e r(n2) r(n3) J ^ ^ ^ 



which is regularized analytically. Therefore 

f(^'='2'=) = . (A.50) 

Other overlap contributions. The remaining overlap contributions in (6.19) are com- 
binations of expansions which already individually lead to scaleless integrals. So it is 
obvious that these contributions are scaleless as well: 



p{h,ic,2c) ^ (ni)ji (n2)j3 (ra3)j3 



-i n J1!J2!J3! 

Jl J2,J3=0 



Bk 



{2pi ■ (2p2 • A;)"2+i2 (A:2)rx3+j3-ii2 

^'^g'^^ V (^l)^l ("2)j. (»3)j34 . 1Y7123^^2y,4 /■ d'^"'^. 

Jl,...J4=0 

00 00 

dk+ {k+y^^^ f dk~ {k-y^^^ 



{-kl)'^3+j3 
— = 0, 

(-A;2 + Qk+ + iO)"i+Ji y (-A;2 + Qk- + iO)^^+h 

—00 

,2.g.7E - (n2),,3 (^3k (-ly^^^ /■ d'^-'fcx(fciy3 



pUcg) ^ /^^^e-^" (^i)ji (n2),,3 (^3k (-1)^"^ /■ 

2z7r'^/2 . ji!---j4! Q"2+i23 _/ 



. . n Q"2+J23 y (_^2 _ ^2)n3+i4 

00 00 

—00 —00 
^(2c,g) ^ Q j^-^ symmetry related to F^^^'^)] , 

,2.ge7E g (ni),, (n2),,3 (^3)^45 (-ly^^M^i'y^ /■ 

jl,...,j5=0 



pih,ic,g) ^ ^'^e-^" (ni),, (n2),,3 (^3),45 (-1?"^^ / d'^-^fc^ 

11 1r,=0 -L' 



00 



^ y (-A?2 + + io)"i+ii y (/c- + io)"2+i3-ii4 - ' 
—00 —00 

^(/r,2c,g) ^ Q [^-^ symmetry related to F^^^'^^'f)] , 

^(ic,2c,,) ^ M^^e^"^ ini)n. (^2),34 (ns),. (-1?^^^ f ^^^^hAgy^ 

2iT,d/2 ,2^ jl!---j5! Q"i2+ii234 y /_p _ 2)n3+i5 

CO 00 

^ y + iO)"l+J2-i35 y (fc- + iO)"2+j4"jl5 ~ ' 
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2iT:d/2 Q«12+il234 



ii,...j6=o 



0, 



(_^2 )n3+i56-i24 J {k+ + iO)"l+J2-J35 _/ [k^ + iO)"2+-'4--'15 
— oo — oo 

p(h,cp) ^ ^(cp) ^ Q jg^j^g expansion] , 

.(lc,cp) _ /^^' e'^"^ (rai)ji (n2)j23 ("3)j45 (-IP^ (m^F" /"H'^-2f: y- 

- 2ivr'^/2 Z_._ ji!---j5! y ^^^^ 

d/fc+ dfe- 



ii,...j5=o 

oo 



~oo 



{k+{k~ + Q) + iO)"'^^' (fc- + iO)"2+J23 (fc+fc- + iO)'^3+J«-J2 
p{2c,cp) ^ Q j^.^ symmetry related to F^^^^^p^] , 
pih,ic,cp) ^ piic,cp) ^ Q [g^^g expansion] , 
pih,2c,cp) ^ p{2c,cp) ^ Q [g^^g expansion] , 

p{ic,2c,cp) ^ /^'^ e^''^ (n2)j34 (n3)j56 (-Ip" f^d~2^ p -^^^ 

2z7r'^/2 ^ ji!---i6! Qni2+ii234 y -Ll±J 

iiv.i6=o 

oo 

r dk+ dk- 

J {k+ + i0)"l+i2-i3 (k- + iO)^2+j4-jl {k+k- + 20)^^3+^56 



~oo 



p{h,ic,2c,cp) ^ p{ic,2c,cp) ^ Q [g^^g expansion] . (A.51) 

These contributions are all well-defined and scaleless through dimensional and analytic 
regularization. 

A. 3. 2 Sudakov form factor: hard contribution 

The hard-region integrals in F^^'^ (6.22), 

DA: 



(A:2 + 2pi ■ A;)"i {k^ + 2p2 • A;)"2 (A:2)"3+j 
are easily evaluated using Feynman parameters, yielding 



(A.52) 



r(ni)r(n2)r(n3+i) 



"1-1 ^n2-l ^na+j-l 
3 



^^""^^^ ^ J dxi dx2 dx3 5(1 - X123) X^i-^ X^^-^ X 





DA; 



2^"l23+i 



((A: + Xipi + X2P2)^ - XiX2<3^) 
„2e e7E g-i7rni23 ( Q2\2-ni23-e-j ^(^123 — 2 + e + j) 

^ ^ ^ r(ni)r(n2)r(n3+j) 

1 

X y dxi dx2 dx3 5(1 - X123) x}-"^^-^-^' . (A.53) 



' V ' 

r(2-n23-e-i) r(2-ni3-e-i) r(n3+i)/r(4-ni23-2e-i) 
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The complete contribution from the hard region is shown in (6.23). In the case ni = n2 
?i3 = 1 this is 



r(l-2e-j) 



j=0 



For the hypergeometric function we use the representations 

and extract the singularity of the t-integral at t — ?• 1: 

1 

-2e 



2Fi(2e, l;l + e;z) = " [^^ " " " 



1 

+ (1 - z)-"^' j dt {1 - ty- 






,2 /'i„2 



= 1 - 2e ln(l -z) + e^ [ln'{l - z) - 2 Li2(z) j + ©(e'^) , (A.56) 

where Li2 is the dilogarithm function. Expanding also the prefactor in (A. 54), the result 
reported in (6.26) is produced. 

A. 3. 3 Sudakov form factor: collinear contributions 

The 1-collinear contribution F'^^^^ in (6.22) reads 

OO / N j 



j=o •' j=o 

D/c 



^7(^2 + 2pi • fc)"i (2p2 ■ A:)"2+i (A;2 - m2)"3-« ' ^^"^''^ 

where ij^y = ((/c^ — m^) + m^)"' in the numerator has been expanded according to the 
binomial theorem. We evaluate this loop integral for general propagator powers ?i-i,2,3 using 
alpha parameters (A. 2): 

y (A;2 + 2pi . A;)"i (2p2 • [k"^ - m2)"3 I JLL T{ni) J ' ' j 



X jBk exp^i(^ai3k^ + 2k ■ {aipi + a2P2)^ ■ (A.58) 
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The loop and alpha-parameter integrations are performed with (A. 4), (A. 2) and (A. 6). 
They yield 

^2e ^-..n,.3 (^2)2-.,3-. (g2)-n. r(n, - u,) Fju,, - 2 + e) r{2 - u, - e) 

r(ni) r(n3)r(2 - n2 - e) 

Plugging this result with n2 — >• n2 + j and ^ — i into (A. 57), we obtain 

^ r(n2_+ J) - "2 - J) / A rjnrs- 2 + e - i) ^^ ^^^ 



j\V{2-n2-e-j) ^ VV r(n3 - i) 



The sum over i is evaluated as follows: 
i 



A t ,\i ^("-13 - 2 + e - i) 



4—n \ / 



r(n3 - i) 



1=0 

1 



r(2-ni-e; 



>!z2 



(-1)^ x-i 



( ^., r(ni3-2 + 6-j)r(2-ni-. + j) 

r(2-ni-e)r(n3) ' ^"^-^^^ 

producing the result which is reported in (6.23). The 2-collinear contribution F^'^'^^ can be 
obtained from F^^'^^ by simply exchanging ni -f-)- n2. 

As explained in section 6, we need analytic regularization for the collinear contribu- 
tions, and we choose to calculate the case ni = 1 + 6, n2 = 1 — S and 723 = 1 in a Laurent 
expansion about 5 = up to the finite order S^. Using the results from above, we obtain 



Q'^\m?) \m?) r(l + (5)r(l -5) ^ V Q"^ 

^ + j) r{26 - j) T{e + 6- j) r{l-e-6 + j) 

j\r{l-e + 5-j) 

1 Y /q2x5 ge7Er(e + <5)r(l-e-5) T{26) 



Q'^\m?) \m? ) T{l-e + 5) r(l 5) 

°° ^m^y (1 - % (6 - 5 ), 

^qy j\{i-26), 



j=0 



2Fi(e-(5,l-5;l-2(5;m2/Q2) 

Using the integral representation (A. 55) the hyper geometric function is expanded about 
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<5 = like in (A. 56): 



2Fi{€-6,l-S;l-26;z) = (1 - z)"' 1 + -5 ln(l - z) 



, fdt 





1-1 + 



tz 



1- z 



+ 0{5^) . (A.63) 



The 1-collinear contribution, expanded in 5, but valid for any e, reads 



1 



+ 



dt 



1 - 1 +t 



Q2 _ 



7E + V(e)-2V'(l-e) 



+ 0{S) , 



(A.64) 



where ipiz) = T'{z)/T{z) is the digamma function. We expand additionally about e = 0, 
keeping the same prefactor {fJ-'^ /Q'^Y as in the hard contribution: 



1 ffi^V 1 1 



2Q2 yQ2 



1 /Q2\ / m2 
+ ln(^)-ln(^l-^ 



+ - In 



2 \ m 
57r2 
12 



+ In 



In 1 



1 2 / 



m 



m 



m 



+ 0{6) + 0{e). 



(A.65) 



Using the symmetry between the two collinear regions upon pi ^ p2 and ni -H- ?i2, the 
2-collinear contribution F^'^^^ can be obtained from (A.65) by replacing 6 — t- —6. 



A. 4 Forward scattering 

A. 4.1 Forward scattering: evaluation with analytic regulators 

In this appendix the evaluation of the contributions (7.19) to the forward-scattering integral 
(7.1) or (7.4) in section 7 is sketched, using the propagator powers ni_2,3,4 as analytic 
regulators. 

Hard contribution. The integral over T^^^I (7.5) is written in a Lorentz-invariant way: 



[h) 



1 f Dk 



1 



+ 



1 



° 2 J (fc2)ni2 \{k'^ - 2pi ■ k)''3 ' (fc2 + 2pi • /c)"3 

1 1 



X 



+ 



(A;2 + 2p2 ■ /c)"4 (A;2 _ 2p2 ■ k^* 



^ + ■,0 -^-;0 ' 
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Bk 



(fc2)"i2 (A;2 _ 2p^ . kyz (A;2 ± 2p2 ■ k) 



n4 



(A.66) 



where terms including the factor {k^ + 2pi ■ have been transformed via k — )■ —k under 

the integral. The two terms F^^^ directly correspond to the diagrams in figure 5 (p. 38). 
They are evaluated in a straightforward way using Feynman parameters and yield 



2e „e7E p-«7i"'T,i234 



(^g2 _ .Q^2-ni234-. 

r(ni234 - 2 + e) r(2 - ni23 - e) r(2 - ni24 - e) 
r(?i3) r(n4) r(4 - ni234 - 2e) 

The complete result reads 

?{^) — ,,2e „e7E p-«7rni234 ^Q2-j2-ni234-e _j_ gi7r(ni234+e) -j 

r(ni234 - 2 + e) r(2 - ni23 - e) r(2 - ni24 - e) 



r(n3) r(n4) r(4 - ni234 - 2e) 



(A.67) 



(A.68) 



Collinear contributions. We also write the integral over T^^'^h (7.7) Lorentz- 
invariant integral using the d-dimensional vector r_|_ = r + [pi — P2)'r'^ IQ"^ satisfying 
Pi,2 ■ = and = — "tJ- With this we express 



r^k + r k^ 
r± ■ k± = —r ■ k -\ = —r± ■ k , 



(A.69) 



using (6.5), and obtain the integral 



DA: 



1 

2 J {k^)^i {{k-r±)^ -2{rl/Q^)p2-k) 
( 1 



"2 



+ 



1 



(A;2 - 2pi ■ A;)"3 (^2 _ ■ k + 2pi ■ k - 2(r\ /Q'^) p2 ■ k) 
1 1 



"3 



+ 



(2p2 • A;)"4 (-2p2 • k) 



724 



(A.70) 



Let us first evaluate only the contribution from the first term in each round bracket 
of (A.70), which is separately well-defined through analytic regularization. Using alpha 
parameters (A. 2), 



1 



2 J (A:2)"i [{k - rx)2 - 2{r\/Q'^)p2 ■ A:)"' (A;2 - 2pi ■ A;)"3 {2p2 ■ ky 



^ n 



^ —iTTrii/2 



\ r(n,) 

X / DA; exp 



da,-a"'-^ e*"2(ri+i0) 



i Oii22,k - 2k 



a2r± + aspi + ( ^ 02 - "4 ) P2 



(A.71) 
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After performing the loop integral with (A. 4) and the a4-integral via (A. 2) we obtain 



1 /U^*^ e'^'^^ g-j7r(ni23+2-e)/2 

2 r(ni)r(n2)r(n3)(Q2 + io)"4 



J dai da2 das a"^ ^ a 



"1-1 ^,"2-1 ^n-i-ni-l ni-2+e 



X exp(i^i^(ri +iO) ) . (A.72) 

V "123 



We can solve this integral by multiplying it with unity in the form of 



oo 



1= fdijdfrj-'^aA, with 0^5C {1,2,3}, (A.73) 

where the sum of alpha parameters in the delta function runs over an arbitrary non-empty 
subset of ai,2,3- Then, changing the order of integration, all alpha parameters are rescaled 
as Oi — )■ rjOi under the 77- integral and the delta function is rewritten as (1/??) ^{^ — YljeS "i)- 
Performing the ?7-integral first using (A. 2) yields 

..2e „e7F, „-j7rni234 / „2 •n\2-ni23-e / ^2 ■n\-n4 ^(72123 — 2 + e) 



00 

X 



2^ V X . V ^ ; r(ni)r(n2)r(n3) 

00 

/ dai da2 das - Yl "i) "2""""' "^""'"^ <2T"^^^' • (^-74) 

Now the three-fold integral is of the form 

da, . . . da„ ^ f 1 - V a,) , "^^"'-f" = J/^^^ " " " , (A.75) 

J^-^ 
valid for any (/> ^ S C {1, . . . ,n}, and we obtain 

^ Ai^^ e^^E e-^^"l234 (^2)2-ni23-e (_g2 _ •Q^-n4 ^^^^ _ 

r(ni23 - 2 + 6) r(2 - ni3 - e) r(2 - n23 - e) , . 
r(ni)r(n2)r(n3)r(4-ni234 - 2e) ^ ' > 

for the contribution from the first term in each factor of (A. 70). Analytic regularization is 
needed for this contribution, otherwise r(n3 — n^) would be ill-defined. 

The contribution from the first x second term in (A. 70), i.e. with (— 2p2 • A:)"""* in- 
stead of (2p2 • A;)"""*, can be calculated from (A. 76) by replacing p2 — )• —p2 and there- 
fore —7- —Q^. (The combination P2/Q'^ remains invariant under this transforma- 
tion.) The remaining two contributions, involving the second term in the first round 
bracket of (A. 70), are obtained from the previous two contributions via the replacement 
— )• r_L + {r'j_/Q'^)p2 — k and ni -f-)- n2 under the loop integral, which leaves the expres- 
sion (A. 76) invariant. The complete result reads 



Fg^l^) = ^2e g.7E g-i7rni234 (^2^2-ni23-. (Q2)-n4 ^ ^i^rn^^ ^ius - n^) 

r(ni23 - 2 + e) r(2 - ma - e) r(2 - nas - e) 



r(ni) r(n2) Tins) r(4 - 711234 - 2e) 



. (A.77) 
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The 2-collinear contribution, 

p{2c) _ 1 f Dfc f 1 1 



2 7 (A;2)"i ((A;-r^)2 + 2(r2/Q2)p^ V(-2pi-A;)"3 (2pi • A;) 

^ +77^^^^ , . , L^ ^ ^ (A-78) 



y{k^ + 2p2 • A:)"4 (A:2 _ 2rx • A: + 2(r2 /Q2) . k - 2p2 ■ k)' 

results from (A. 70) by exchanging pi -H- —p2 and 713 -H- 714. It is thus obtained from (A. 77) 
via 713 -f-)- 714 and reads 



FfS^'^) = ^2. g.7E g-i^ni234 (^2)2-ni24-e (Q2)-n3 ^ e^^"^) ^(714 - ng) 

r(rii24 - 2 + e) r(2 - tiu - e) r(2 - 7134 - e) 



^0 



T{ni) T{n2) T{ni) r(4 - 711234 - 2e) 

I 



(A.79) 



Glauber contribution. The integral over T^^h (7.11) reads 



2 2i7r'^/2 7 (_^2)ni (_(^^_^-^)2)"2 



00 



X / dk+ \ \ + 



-00 
00 



{-Qk+ - k{ + z0)"3 iQk+ -k{ + 2r± ■ k± + iO) 



dk' I 1 + I . (A.80) 

\{Qk~ -kl + iO)^^ {-Qk~ - kl + 2r^ ■ k^ + iO)^^ I 

Analytic regularization makes the /c^-integrations over the individual terms well-defined 
and allows to assume Re 713^4 > 1. Then we can close the integration contours at ioo or 
—zoo, choosing separately for each term the side where no pole and no branch cut lies within 
the closed contour. So all k~^- and A;~-integrals vanish as scaleless integrals, cf. (A. 35). The 
fej_-integral is regularized dimensionally. Therefore, 

= . (A.81) 
Collinear-plane contribution. Integrating over r(S'^^/ (7.13) yields 

_ 1 /■ ,_2^ f dk+dk 



2 2i7r'^/2 J ^ J {k+k- + i0)"i2 



1 1 



+ 



\^ {k+ {k- -Q) + iO) {k+ {k~ + Q) + iO) J 

(y{{k+ + Q)k- + iO^' ^ {{k+ - Q)k~ + iOy") ■ ^^'^^^ 

As the integrand is independent of k± , the A:_L-integration is scaleless through dimensional 
regularization, cf. (A. 38). The fc^-integrals are regularized analytically. So 

F(S'^)=0. (A.83) 
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Hard— collinear overlap contributions. The hard-l-collinear overlap contribution 
reads 



I r Bk / 1 



° 2 J (A;2)"i2 V^(fc2 _ . A;)"3 ' (A;2 + 2pi • fc)"3 

1 



It corresponds to Fq^'^^ (A. 70) with r_|_ = 0. According to (A. 72), we get 

' - ' r(,.)r(,.)r(„3) (<« + + <-<^ + 

ni — 1 712 — 1 n3— n4 — 1 n4— 2+e 



X y dai da2 das a"^ ""^ ^ a 



'123 



,,2e e7E p-«7r(ni23+2-e)/2 VO — ri^^A—e} (' 

(Q2)n4 V ^ r(n3)r(2-n4-e) 7 ^ 



= 0. (A.85) 

This integral is scaleless by dimensional regularization. Analytic regularization is only 
required if the aa-integration is performed first. Similarly, via the symmetry pi -H- —p2 
and n3 -H- 77-4, 

Fq^'*'^'^) = . (A.86) 

Collinear collinear overlap contribution. The 1-collinear 2-collinear overlap contri- 
bution reads 



° 2 J (A;2)ni ((^ _ r^)2)"2 \^(_2pi • fc)"3 ^ (2pi • k^s 



^ .w. )• (A.87) 



X2p2-kYi {-2p2-kY^ 
The first term of each round bracket yields the contribution 

X y DA; exp^i[ai2A;^ — 2A; • (02^^ + a^pi — 04^2) + ^O] 

00 

^ 1 ^2.g.7Ee-^^("l23+2-e)/2 . 

2r(ni)r(n2)r(n3)(Q2 + iO)"4y % 



00 

X / dai da2 q"^"'^ Q2^~^ Q;i2~^"'"'^ 6xp ( i ^ ^ (^1 + ^0) ) • (Ai 

J V "12 y 
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The ai^2-integrations are regularized dimensionally, while the as-integral is scaleless by 
analytic regularization. The remaining contributions from (A. 87) are obtained via pi — >• 
—pi and/or p2 — )■ —p2 and vanish as well. Therefore, 



Other overlap contributions. All other overlap contributions in (7.19) are expansions 
of scaleless integrals evaluated above, so they are scaleless as well: 



1 f Dk f 1 1 



^0 



^0 



2 } (A;2)'^i2 \{-2px-kY'i (2pi • A;)"3 

1 +__L_i=„, 



Xlp2-kY'^ {-2p2-kY^ 
1 [x^^ e^^E r d'^-^A?^ 



2 2m'^/2 ] {-kl)ni2 

X / dfc+ ' ^ 



-oo 
oo 




X / dA;- + =0, 

J \{Qk- -k] +iO)"4 ' - ' 

—oo ^ 

ilc,g) _ 1 f^^' e'^^ 1 f d^-^ki 



oo / 

X / dk+ ( 1.- + ^ 



(-(5/c+ - A;^ + i0)"3 (Q/c+ - /c^ + 2fj_ • kj_ + i0)"3 

-oo ^ 

oo 

X / dA;~ ( — ^ h — ^ — 1=0. 

\(A;--FiO)"4 io)"4 



Fq = [via symmetry related to F^ 



1 /i^^e^^E 1 y d-^-a^^ 



00 / 

X / dk+ ( + ^ 



{-Qk+ - A;2 i0)"3 (Q/c+ - A:2 iO)"^ 
—00 ~ • 
00 

X / dA;~ I r h ^ — 1=0, 

—00 

p[h,2c,g) _ Q j^.^ symmetry related to Fq^'^^'^^] , 
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p(lc,2c,g) i 



oo 

X fdk+( , , , ^ + ^ 



(-A;+ + iO)"3 (fc+ + io)"3 



-oo 

oo 



X f dk~ ( 3— h — ^1=0, 

— oo 

_ oo 



— oo 

oo 



X f dk (- — ^ \- —- — ^ — 1=0, 

— oo 

^ih,cp) ^ p(cp) ^ Q jg^^^ expansion] , 

oo 

^(iccp) ^ 1 /^^^e-^^ J_ /■ ,_2^ /■ d/c+dfc- 
° 2 2i7r'^/2 Q"4 y -L y (fc+yt- + iO)"i 

— oo 

/ 1 1 

+ 



(A;+(A;- -Q) + iO)"' (A:+(A:- + Q) + iO 

+ TTT^ =0. 



(A;-+iO)"4 (_A;- + io)"4 
^^(2c,cp) _ Q j^.^ symmetry related to i^^j^'^''^^^ 
p(h,ic,cp) _ p(ic,cp) _ Q [game expansion] , 
p{h,2c,cp) _ p{2c,cp) _ Q [same expansion] , 

oo 

piic,2c,cp) ^ 1 fi'-e^^^ 1 /-^.-s^ /■ dA;+ d^ 



° 2 2i7r'^/2 Q"34 y y +^o)"i2 



-oo 



1 1 \ / 1 1 , 



^^(h,ic,2c,cp) ^ p{ic,2c,cp) ^ Q jg^^^ expansion] . (A.90) 
All these integrals are well-defined via dimensional and analytic regularization. 



A. 4.2 Forward scattering: evaluation without analytic regulators 

Here the evaluation of the contributions (7.19) to the forward-scattering integral in section 7 
is repeated with analytic regularization switched off. This means that the propagator 
powers are set to the fixed values ni = n2 = ^3 = 714 = 1 from the start. 
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Hard contribution. The evaluation of the hard contribution in appendix A. 4.1 is also 
valid without analytic regularization. From (A. 68) we get 

"(W wJ ^ ^ n^) • ^ 

CoUinear contributions. The 1-collinear contribution without analytic regulators reads 

oo 



2Q 2!ir<i/2 J J lk+k~ - e + lO) - - (fcx - rl)2 + tO) 

— CJO ^ \ u / 

1 1 \ 

+ 



k+ik~ -Q)-kl + iO k+{k~ + Q) - r+k~ -kl + 2r^ -kx+iOj 



k- + iO -k-+iO 

The factor in the last line of (A. 92) vanishes everywhere under the ^"-integral except for 
the pole at = which lies below the real axis in the first term and above the real axis 
in the second term. The only remaining contribution arises from the closed integration 
contour around this pole, i.e. from the residue at k" = 0, which allows us to express 

^ + —r^ — - = 1—^—77 - 1— ^— 7T = -2i7r 6{k-) . (A.93) 



k- + iO -k~ + iO k~ + iO k~ - iO 
The 1-collinear contribution becomes 



^(1,) 1 /i^^e^^E f 
rn ~ 



2Q vri-^ J p {k^ - fx)2 



dk+ I ^ + ^ ^ |. (A.94) 

V -Qk+ -kl + iO Qk+ -kl + 2ri_ ■k± + iOl 

The A;~''-integral converges at |A;+| —7- 00 because the leading behaviour for large \k~^\ is 
cancelled between the two terms, such that the integrand falls off like l/{k~^)'^ towards 
infinity. We close the integration contour at +ioo or at —ioo, taking into account the 
residue of the one pole within the contour. In general, we get 

00 

, /I 1 \ _ 2i7r 

[-Qk+ + A + iO^ Qk+ + B + io) ~~~Q' ^ ' 

—00 

for real A, B and Q > 0, which is in particular independent of A and B. The fc^-integral 
is solved using a Feynman parameter. 



k\{ki_-f^Y J J (^(k± - xf±)^ + x{l - x)fl) 



2 
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1 

r(l + ^) /d. (x(l - .) ri)-'-' = ^ ^^ll^^l'^ ■ (A.96) 



The complete 1-collinear contribution reads 



^7^ 



ee7Er(l + e)r2(-e 



r2Q2 \^f2 J T{-2e) 
The 2-coninear contribution, 



2Q 2i7r'^/2 7 J ^k+k- -kl + iO){k+{k- -ro)-ik±-r±y+iO) 

1 ^ \ 

+ Q)A;- -kl + iO (fc+ - Q)A:- - k+r^ -kl + 2f^ -k^+io) 



-A;++zO k++iO 

is obtained from -Fg"'^'^^ (A. 92) by replacing k~^ —k~ and r^" — t- — = '"o'- So the same 
result is reproduced, 

^(20 _ ^(10 _^(llY e-^-r(l + e)r^(-e) 
^0 " ^0 " rl Q2 ^^-2 j r(-2e) ' ^^'^^^ 

Glauber contribution. The Glauber-region integral (A. 80) for ni = n2 = = n4 = 1 

reads 

^(3) 1 n^'e'^^ f d'^-^k 
^0 ~ - - 



2 2i7r'^/2 J p {k^_ - fx)2 



00 / 

X I dk+ I ^ + ^ 



-00 
00 



-Qk+ -kl + iO Qk+ -k']+ 2r± ■ k± + iO 



X / dA;- ( ^- + — |. (A.lOO) 

J V Qk~ -k] + iO -Qk- - fc^ + 2f^ ■k±+iO 

—00 ^ ' 

The fc^-integrals converge because the leading behaviour for \k^\ — t- 00 is cancelled between 
the two terms of each round bracket. All integrals needed for (A.lOO) are known from (A. 95) 
and (A.96). The result is 

-^l^J — n^) — • ^^-'"'^ 

Collinear— coUinear overlap contribution. 

00 

,(ic,2c)_ 1 ^^'e^^E f [■ dk+dk~ 



2Q2 2i7r'^/2 J J (fc+fc- -p + iO) (k+k- - (k^ - f^y + iO) 

—00 



-A;+ + iO k+ + iOj \k-+iO -k-+iO 
■' ^ 



-2't7r(5(fc+) -2't7r(5(fc-) 
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where the two round brackets are replaced according to (A. 93). Using (A. 96), the 
l-colhnear-2-coninear overlap contribution yields 



-^UJ — — ■ ^ ^ ^ 

Collinear Glauber overlap contributions. The 1-collinear-Glauber overlap contri- 
bution 

_ oo 



2Q 2i7r'^/2 J k\ {k± -r^Y J \k~ + iQ -k~ + iO 




^ + ^, , ro r I ' (A- 104) 



-Qk+ -k] + iO Qk+ - fe^ + 2fj_ • A;_l + iO 



its symmetric counterpart F^^^'^^ (obtained via A;"*" -H- —k ), and the threefold overlap 
contribution 



oo 



-oo 



oo 



X I dk- { - — h —r^ 1 (A.105) 

' ^ k- +iO -k- +iO ' ^ ' 



-oo 



are evaluated using the relations (A. 93), (A. 95) and (A. 96). They yield identical results, 

p(lc,g) _ {2c,g) _ p{lc,2c,g) _ ( f^^V e'^I^T{l + e) T^{-e) /A Iflfi^ 

-WqA'^J — — ■ ^ ^ 

Contributions involving the coUinear-plane region. The contributions F^^^^ (A. 82) 
and F^''^'^\ p(Kic,cp) ^ j,(h,2c,cp) ^ p{ic,2c,cp) ^ j,{hM,2c,cp) ^^ ^^^ .^^^^^^ 

the scaleless integral / d'^~'^k± (with no further integrand) which vanishes through dimen- 
sional regularization. This is still true without analytic regularization, when the propagator 
powers are fixed to ni = n2 = = = 1. But the /c^-integrals in these contributions 
are problematic without analytic regulators. The integrands fall off sufficiently fast at 

— )• oo, but the collinear-plane expansion makes poles above and below the real /c^-axis 
collapse at k^ = 0. So the integration contour is pinched between poles which are infinites- 
imally close to each other. These /c^-integrals are regularized analytically in the evaluation 
of appendix A.4.1, but here they are singular. 

We may regularize the pinching of poles in the fc^-integrals without resorting to an- 
alytic regularization by keeping a finite (instead of infinitesimal) imaginary part iS (with 
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(5 > 0) in all propagators, e.g. 



oo 



2 2i7r'^/2 j ^ J {k+k- + iSy 



1 1 

+ 



k+{k~ -Q) + i5 k+{k- +Q)+i5 

^ + ..^ , ... ) • (A.107) 



{k++Q)k~+i6 {k+ - Q)k~ + i5 ) ' 

When 5 is fixed at some finite positive value, the poles above and below the integration 
contours are separated from each other and the collinear-plane integral is well-defined. 
The result from the /c^-integrations is singular in the limit (5 — t- 0, but it is multiplied 
by the scaleless /i;j_-integral, so the complete contribution vanishes. This is true for all 
contributions involving the collinear-plane region: 

_ rp{h,cp) _ p(lc,cp) _ p{2c,cp) _ p{h,lc,cp) _ p(h,2c,cp) _ p(lc,2c,cp) 
^pih,lc,2c,cp) (A. 108) 

For consistency, the same regularization by a finite imaginary part i5 should be applied to 
the contributions from the other regions as well, but there the limit (5—7-0 (with 5 > 0) is 
regular when it is performed before the limit d — t- 4 from dimensional regularization. 

In addition, as we have already noted after the identity (7.19), all contributions involv- 
ing the collinear-plane expansion T^'^^^ cancel each other in pairs at the integrand level, 
according to (6.21), because the hard expansion T^^^ does not change the integrand when 
it is applied in addition to T^'^p^: 

^^^(cp) _ p.(h,cp)-^ _ ^-p,{lc,cp) _ p,{h,lc,cp)-^ _ ^■p,{2c,cp) _ p{h,2c,cp)^ 

^ ^p{lc,2c,cp) _ p(h,lc,2c,cp)^ ^ Q ^ ^^^^Qg^ 

So even if individual collinear-plane contributions were singular, their sum would vanish at 
the integrand level. Also, the collinear-plane contribution Fq^^^ is parametrically of order 
(^1)^ {Q^)~^ ('''•21) and thus suppressed by two powers of f]_/Q'^ with respect to Fq^^\ 
F^j^^^ and Fq^^ . So it does not contribute to the leading-order result Fq . 

Overlap contributions involving the hard region. The overlap contributions involv- 
ing the hard expansion and any of the collinear or Glauber regions read 

oo 

p{Mc) ^^^l^[.d^2r f dfe+dfc- / 1 1 

^0 = o._rf/2 / d / 77T7Z — , „-n + 



2Q 2i^'^/2 J ^k+k- -kl+iO)^ \k~ +iO -k- + iO 



1 1 

+ 



k+{k- - Q) - kl + iO k+{k- + Q) - kl + iO 



-76- 



oo 



1 ^u^^e^TE fd'^-^k± /•,,+ / 1 1 



i^o^''''^ = =: ^^-TTTT / ^^^r^ / dk+ I —^^ + 



2 2i^'^/2 7 (^^2)2 7 \Qk+-kl+iO Qk+-kl + iO 



X / dA:- I Ir- + ^ 



Qk- - fc2 +iO -Qk- -k^. +iO 



oo 



2Q2 2ii<i/2 J {k*k--kl+iO) 



-OO 



1 , 1 ' 



A;+ + iO A;+ + iO / V fc" + «0 -k~ + iO 



^{/i,lc,g) ^ ^ ^- - ; " / ^f^+ I : , : 

° 2Q 2i^'^/2 7 (^2)2 J y_Qk+-kl + iO Qk+-kl + iO^ 



1 t?^e^"i^ fd'^-^k± 



oo 

X I dk- { —3—- + ^ 



k~ + zO -k- + iO 

_ oo 



2g2 2m'^/2 J (^2)2 7 \^_/s++io fc++iO 



oo 



X / dfc- ( - — i + — -^^ ) , (A.llO) 



omitting the symmetric integrals F^''^^^ and p^''^'^'^^ obtained from the integrals above via 
(Ic) — )• (2c) with /c"^ -H- — /c". The A;^-integrals are evaluated using (A. 93) and (A. 95). All 
contributions in (A.llO) identically yield 

p(/i,lc) _ Tp{h,2c) _ rp(h,g) _ p(h,lc,2c) _ p{h,lc,g) _ j^{h,2c,g) _ p{h,lc,2c,g) 

= ^ li^- / "±^±1. = , (A.lll) 
vri-^ 7 (fc2)2 ^ ^ 

which is a scaleless integral through dimensional regularization. 



B Expansion by regions with a finite boundary 

This appendix illustrates the behaviour of the expansion by regions and the formalism 
of sections 3 and 5 when a finite integration boundary is involved. Consider the one- 
dimensional integral 

B 



F = dkl with / = — and < m < A/ < S . (B.l) 

' {k + m)"' {k + M)l^ ^ ' 
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B.l Direct evaluation with a finite boundary 

The integration domain D = [0, B] in (B.l) involves tlie finite, non-zero boundary B, which 
introduces an additional scale into the problem. Let us define two regions, 

• the soft region (s), characterized by A; ~ m, with the expansion 

ii=0 ■'^ V -r ; 

converging absolutely for A; € = [0, A] with m <C A ^ M, 

• and the hard region {h), characterized by A: ~ M, with the expansion 

~i^o ' ' -^2' A:"+J2 + ' ^ ^ ^ 

converging absolutely for k G = (A, S]. 
We have Dg D Dh = i/i, U Dh = D and 

oo 

ilj2 = 

71.79 „• I „• I ji/ffl-i-ii lo — 1 1 ■ V • y 



Also the integrals 





B 

(h) _ {a)j2 {-my^ f dk 



F 



j^l J k°'+n {k + MY ' 





are well-defined if we use the denominator power a as analytic regulator, choosing its value 
for each of the integrals F-'^^ and F^^'^^ such that they converge at A; — >■ 0. So the conditions 
1, 2 and 3 of the formalism in section 3 hold. But we have a problem with the convergence 
of the summations as required in condition 4: The series X^°^=o ^'fx ^'^d X]j^=o ^'f/jl (^°^ 
fixed 22) diverge, because these two integrals involve contributions from k ^ B producing 
powers {B/Uy^ (with B/M > l).*' 

Note, however, that the series ji) are convergent, although 

contributions (m/ky^ with — >■ are involved. But the analytic regularization prevents contributions from 
the zero-boundary fc = and makes the integrals yield results scaling only with fc ~ M and k ^ B. 
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With F^*) and F^'^'^') being divergent and therefore condition 4 violated, we cannot 
naively claim the master identity (3.22), here F = — _p(*.'^)^ to be true. However, 

the difference F^^') — F^^'^^ is finite when the terms are combined before summing over ji. 
This can be seen by writing 



_B oo oo 

J {k + m)° ~ J {k + m)" J {k + mY 

B 



oo 



ji! r(a - 1 - ji) ^ (a),, i-my^ f dk 



using the power a as analytic regulator also at /c — ?• 00 and expanding the second term 
safely for m <^ B < k. The two contributions combined yield 

p(s) _ pis,h) ^ i/3)n r(a - 1 - ii) /_^yi 
i2=o ^ 

^ / B 00 \ 

(«).2(-^^yM /• dk f dk \ 

7i!M/3+ii ^ 72! \J k°+h-h ^ J A:"+i2-ii I • ^ ' ^ 

J2=o \o B J 

The sum of the two integrals in the second line of (B.7) is zero (using analytic regulariza- 
tion), because 

^ dk 5i-"+ii-J2 , . . 

ua+io-h = 1 r~ ~ for Rea < 1 +ji - j2 , 

/■ dA- Rl-a+ii-i2 
T:^Tn^ = -r, — for Rea>l+ii-i2, (B.8) 

B 

or simply because the two integrals add up to a scaleless integral from to 00. 

The upper limit of the integrals Fj^^ in (B.5) can be extended to 00 like in (B.6), and 
the second term is evaluated in analogy to (B.8). In total the result is 





00 



00 \ 00 

jl=0 \ J2=0 / j2=0 



^MP a-l ^ \m) (2-a)ji 



ii=o 



1 r(l-a) ^ /rriy2 r(a + /3-l+i2) 



+ M^T^ r(/3) \mJ ja! 

1 ^ / MVV my2 (/3),, (g),. 



I V bJ ji!j2!(a + /3-l+ii2 
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The series expansions in (B.9) converge and reproduce the exact result of the integral F 
in (B.l). Note that due to the finite integration boundary oXk = B the individual integrals 
F^^ and F^^^ in (B.5) do not yield homogeneous functions of m, M and B, and therefore 

the overlap contributions F- '• are not scaleless. Additional expansions of Fi^^ and Fi''^ as 
in (B.6) are needed to arrive at the form (B.9) where all summation terms are homogeneous 
functions of m, M and B. 

Individual terms in the result (B.9) are singular when a is a positive integer or when 
(a + /3 — 1) is zero or a negative integer. One can check, however, that all these singularities 
cancel between the terms in (B.9) such that the total result is finite as it should be. 

The difficulties with the convergence of the expansions F^'*) and F(^''^) may be avoided 
if one is only interested in a leading-order expansion. According to (3.27), the leading-order 
approximation to F is given by 

B B B 

M j_ Uh) M) 1 f dk f dk I fdk 



Fo-Fo +Fo Fo^^ - MP j {k + mr^ j k-{k + MY J k- ' ^^'^^^ 



with the zeroth-order terms from (B.5). It does not require the convergence of the full series. 
The first two contributions to the leading-order approximation Fq are not homogeneous 
functions of m, M and B; they may be further approximated using (B.6) and keeping 
only the leading term of the additional expansion. This leads to a different leading-order 
approximation with homogeneous contributions: 

oo oo oo oo 

1 f dk f dk f dk 1 fdk 

+ 



M/3 j {k + mY J k'^{k + MY J k°'+l^ Ml^ J k» 

BO 

1 1 1 r(l-a)r(a + /j-l) 1 1 



which corresponds to the leading-order term (with ji = j2 = 0) in (B.9). 
B.2 Evaluation by extending to an infinite boundary 

There is an alternative way to treat the integral (B.l) with finite boundaries in the ex- 
pansion by regions, which, however, requires the formalism for non-commuting expansions 
introduced in section 5. 

Let us write the original integral as 

oo 

dki with /=- (B.12) 

(A; + m)" (fc + M)/5 ' ^ ^ 



with infinite upper boundary of the integration domain D = [0,oo), expressing the restric- 
tion k < B through the Heaviside step function 6. With a third factor and the new scale B 
in the integrand /, we add a third region. Let us keep Ds = [0, A], but redefine = (A, A] 
with m<cA<CM^A<Ci?. So, in the soft and hard regions, we have k < A <^ B, which 
allows us to use 

e{B-k) = l, k<B, (B.13) 
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there. ^ So the soft and hard expansions of the new integrand, T^'*) and T^^\ as well as 
their combined double expansion T^*''*) exactly correspond to the known expressions (B.2), 
(B.3) and (B.4). Additionally we now have 

• the ultrahard region (uh), characterized by A; ~ i?, with the expansion 
ji,h=o 

converging absolutely for k G = (A,oo). 

We have chosen the new convergence domains such that DsCiDh = DsCiDuh = DhCiDuh = 
and Dg U U D^^ = D. The ultrahard expansion commutes with the hard expansion, 

oo 

rp{h)j^{uh) _ rp{uh)j^(h) _ rjn{h,uh) _ rp{h,uh) 

ii j2=o 

but not with the soft expansion: 

rp{uh)rp{s) — rp{uh-(—s) Ji(s,h) 

rp(s)rp{uh) _ rp{s<^uh) _ rp{h,uh) -j^g-^ 

Therefore we need the formalism for non-commuting expansions as developed in section 5. 
The subset of regions with non-commuting expansions is i?nc = {s,uh}, and condition 2a 
(p. 25) holds with Rc = {h}. Also condition 5 (p. 28) is fulfilled because the hard expansion 
does not further change the integrands which are already doubly expanded with T^^^ and 

'-p{uh) . 

rp[uh<—S,h) rp{uh^s) rp(s<—uh,h) rp{s<—uh) -j^'^j 

The integrals needed within this formalism read 

oo 

pis) ^ W)n (-ly^ / <ikkn ^ 1 / m yi _J_ {(3),, 
h j^lMl^+n J {k + mY m°'-'^M'^\MJ a-l{2-a)j^' 



oo 

p{h) _ {a)j2 (-m)J2 f dk 



32 



J2! J k^+n {k + MY 


1 /my2 r(l -a) r(a + /3- 1 +^2) 



M^+Z'-i \M) r(/3) 32 



h,32 jl\j2^- J fc"+/3+ii2 







1 f My^/_my2 mjA»)j2 



V BJ \ B) ji!j2!(a + /3-l+ji2) 



■^Formally this is the expansion e{B - fc) = e{B) - k5{B) + \ 5'{B) + ... = 1 with 61(B) = 1 and 
5^^\B) = for B > to all orders j. 
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oo 

hJ2 ji\i2^.Ml^+^'^ J k'^+h-ji ' 





oo 



^h,n -f-! J k^+P+n2 (^-1^) 



i.e. all combinations of expansions which commute with each other. All integrals in (B.18) 
are well-defined via the analytic regulator a (condition 3) , and the summations over ji , j2 
converge absolutely for m <^ M <^ B (condition 4). 

According to (5.14), the integral can be expressed as 

p — p{s) _|_ pih) _|_ p{uh) _ pis,h) _ p{h,uh) 

j-,{uh-(—s) j-,{uh-(—s,h) Yp{s-(—uh) Yp{s-(—uh,h) /-p q\ 

~^[uh] '^^[uh] ~ [s] W ■ [O.iy) 

As condition 5 holds, the extra terms in the second line of (B.19) cancel with each other. 
So according to (5.18), 

p — pis) _|_ pih) _|_ piuh) _ pis,h) _ p{h,uh) _ (B.20) 

The first three terms with single expansions are series of homogeneous functions of m, 
M and B, so the overlap contributions in the last two terms yield scaleless integrals and 
vanish, see (B.18). The final result shortens to 

F = + F^^^ + , (B.21) 

reproducing the result (B.9) from the previous calculation, when the expressions in (B.18) 
are added together. 

Let us have a short look at the extra terms in (B.19) which we have dropped: 

oo 

^{uh^s) ^ piuh^s,h) ^ {f^)n (Q)j2 (-IP' (-m)J2 r dk 

j2,jlluh] j2,jiluh] ji! M^+Ji J 

A 

1 ( Ay^/_my2 (/3),,(a),, 



M) V A/ ii!j2! (1 - a +Ji - J2) 



A 



^(^s^uh) ^ pis^uh,h) _ (a)j2 {-My^ {-my^ f dk 



J1J2H ji,j2ls] ii!j2! J A;°+^+-''i2 



1 f MY f my 2 



- 



^a+p-iy V XJ ji!j2!(a + /3-l+ji2) 



(B.22) 



The summations of these terms over ji diverge, because A/M ^ 1 and M/X ^ 1. So 
by choosing to expand the integral F (B.l) via extending its upper boundary to infinity 
and introducing a third region, we still cannot completely avoid the problem of diverging 
series expansions encountered in section B.l. But at least we have isolated the divergence 
problems in terms which disappear from the final result (B.21) in a well-understood way. 
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Note that also other series expansions which appear in intermediate steps of the deriva- 
tion, in particular V • F^^''^\,, and V ■ F^^'^'ll, , diverge because the domain Dh has bound- 
aries touching both the soft and the ultrahard domain. So we can only sum the expansions 
over ji up to oo once the terms have been combined into the integrals (B.18) performed 
over the complete integration domain. For those the convergence of the series expansions 
is required by condition 4 in section 3. 
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